Competing order and nature of the pairing state in the iron pnictides 
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We show that the competition between magnetism and superconductivity can be used to de- 
termine the pairing state in the iron arsenides. To this end we demonstrate that the itinerant 
antiferromagnetic phase (AFM) and the unconventional s^ sign-changing superconducting state 
(SC) are near the borderline of microscopic coexistence and macroscopic phase separation, explain- 
ing the experimentally observed competition of both ordered states. In contrast, conventional s^^ 
pairing is not able to coexist with magnetism. Expanding the microscopic free energy of the system 
with competing orders around the multicritical point, we find that static magnetism plays the role of 
an intrinsic interband Josephson coupling, making the phase diagram sensitive to the symmetry of 
the Cooper pair wavefunction. We relate this result to the quasiparticle excitation spectrum and to 
the emergent S0(5) symmetry of systems with particle-hole symmetry. Our results rely on the as- 
sumption that the same electrons that form the ordered moment contribute to the superconducting 
condensate and that the system is close to particle-hole symmetry. We also compare the suppression 
of SC in different regions of the FeAs phase diagram, showing that while in the underdoped side 
it is due to the competition with AFM, in the overdoped side it is related to the disappearance of 
pockets from the Fermi surface. 



I. INTRODUCTION 

The recent discovery of iron arsenide 
superconductors^"— brought renewed interest in the 
research of high-temperature superconductors. With 
transition temperatures Tc of more than 50K in some 
cases, these compounds present a very rich phase 
diagram displaying superconducting (SC), antiferro- 
magnetic (AFM) and structural order—"—. In some 
pnictides, such as LaFeAs (Oi_xFx), PrFeAs (Oi_xFx), 
(Sri_xNax)Fe2As2, and (Bai_xKx) Fe2As2, the com- 
peting SC and AFM phases seem to be separated 
by a first-order transition and can only coexist in 
phase-separated macroscopic regions of the sample^i^rJ^. 
However, in other compounds, like Ba(Fei_xCox)2 As2 
and possibly^ SmFeAs (Oi_xFx), local probesi^"— as 
well as bulk measurementsi^ii^"— demonstrated that SC 
and AFM coexist homogeneously. This coexistence, how- 
ever, is characterized by a competition of the two ordered 
states: neutron diffraction experiment a^^i^^ revealed the 
dramatic suppression of the magnetization below Tc, 
to the extent that reentrance of the non-magnetically 
ordered phase sets in at low temperatures^Si. 

Experiments have also demonstrated the itinerant 
character of the magnetically ordered phase in the pnic- 
tides. In particular, optical conductivity measurements 
show a considerable Drude weight as well as a pronounced 
mid-infrared peak below the Neel transition temperature 
T/v, consistent with the itinerant pictureSii?-^. Further- 
more, band structure calculations reveal that the crys- 
talline field is unable to significantly split the energy lev- 
els in order to localize 3d electrons^^. Also, several theo- 
retical models demonstrate the adequacy of the itinerant 
description^ir— . Therefore, in the iron arsenides, the 
same electrons that form the superconducting conden- 
sate seem to be the ones that contribute to the ordered 
moment. 



The interplay between AFM and SC has been 
investigated in many contexts^Sr^, including the 
pnictides^-!^^"— . In this paper, following results from 
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we mves- 



our previous work=ii as well as from Refs 
tigate in detail the connection between the competition 
of these two phases and the pairing symmetry of the SC 
state. We demonstrate that coexistence between SC and 
itinerant AFM in the pnictides for temperatures close to 
Tjv — Tc is only possible if the pairing state is uncon- 
ventional, as proposed by models with purely electronic 
pairing mechanisms^"-^. In particular, using a mean- 
field Hamiltonian for the competition between AFM and 
SC and expanding the microscopic free energy in powers 
of the order parameters, we show that a conventional s++ 
SC state does not allow a coexistence regime to be estab- 
lished around the point where the Tn and Tc lines meet, 
even for extreme values of the band structure param- 
eters. Meanwhile, the unconventional s^ state, whose 
gap function changes sign from one Fermi surface sheet to 
the other, may or may not coexist with AFM, depending 
on the details of the band structure dispersion relations. 
Specifically, for the parameters of Ba (Fei_xCox)2 As2, 
we find that the s^ state coexists with itinerant AFM, 
while the s++ state does not (see figure [Ij. 



Our results rely solely on the general assumptions that 
the magnetism is itinerant and that the band structure of 
the iron pnictides is not far from particle-hole symmetry, 
consisting of two distinct sets of Fermi surface sheets^: 
hole pockets located at the center of the Brillouin zone 
and electron pockets displaced from the zone center by 
the magnetic ordering vector Q. Additional details of the 
band structure, the dimensionality of the system or the 
presence of intraband pairing interactions do not change 
the conclusions. On the other hand, for a localized AFM 
state, the free energy expansion reveals that coexistence 
is easily attained, which is difficult to reconcile with the 
observation of phase separation in some compounds. Our 
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Figure 1: Phase diagrams of Ba(Fei_xCox)2 As2 for a su- 
perconducting s^~ state (a) and an s"*"^ state (b), obtained 
by numerically solving the gap equations. The green region 
denotes homogeneous, microscopic coexistence, whereas the 
dark red region denotes heterogeneous, macroscopic coexis- 
tence. The band structure parameters are discussed in Sec- 
tion IV-B. 



analysis also indicates that the onset of SC has little effect 
on localized moments, which is at odds with experimen- 
tal observations as well, giving further evidence for the 
itinerant magnetism in the pnictides. 

We also investigate in detail the origin of the strong 
dependence of the phase diagram on the symmetry of 
the Cooper pair wavefunction. Expressing the Ginzburg- 
Landau coefficients in terms of Feynman diagrams, it be- 
comes clear that the static staggered magnetic moment 
m plays the role of an intrinsic interband Josephson cou- 
pling. Specifically, it corresponds to a term in the free 
energy of the form: 



Ej (xm^lAillAalcos^ 



(1) 



where 9 is the relative phase between the SC order pa- 
rameters of the two Fermi surface sheets, Ai and A2. 
Thus, the coexistence state in some iron arsenides nat- 
urally carries information about the relative phase of 
the Cooper pair wave functions, which are usually ac- 



cessible only through intricate and delicate interference 
experiments^"—. 

The quasiparticle excitation spectrum is substantially 
different for distinct pairing symmetries. For the spe- 
cial case of particle-hole symmetric bands, the system 
with competing magnetism and 5"*"+ pairing have two dis- 
tinct positive eigenvalues i?k with E^ — ^l^ + {M ± A) , 
whereas for s^ pairing the positive eigenvalues are de- 
generate: ^k = ^k + ^^^ + A^. We find that this 
special form of the excitation energy in the s^ case 
implies that all quartic (and higher order) Ginzburg- 
Landau terms must depend solely on the combination 
(Af^ -I- A^), which is the root of the S0(5) symmetry of 
the system and, ultimately, what leads to the conclusion 
that AFM and SC are in the borderline of coexistence 
and mutual exclusion. Inclusion of local moments add 
new terms to the free energy, removing the system from 
this borderline regime. 

In this paper we also investigate some specific prop- 
erties of the coexistence AFM-SC phase. We find that 
the effects of Coulomb repulsion on the magnetic super- 
conducting phase are basically the same as in the case 
of a pure multiband superconductor. Particularly, the 
s^ state is remarkably stable with respect to an uni- 
form Coulomb repulsion. We also studied analytically 
the shape of the reentrant Neel transition line inside the 
SC state for low temperatures. At the mean-field level, 
the finite SC gap introduces an overall energy scale that 
causes the Neel line to have a divergent slope as T — > 0. 
Quantum fluctuations, which are relevant only in a very 
small region, end up suppressing the reentrant behavior. 

The competition between AFM and SC explains the 
suppression of Tc in the underdoped side of the FeAs 
phase diagram. Here, we investigate the suppression of 
SC on the overdoped side as well. We find that the 
changes in the Fermi surface with doping are crucial to 
kill SC, in agreement to ARPES measurements^"— . In 
particular, the vanishing of one of the Fermi surface pock- 
ets marks the onset of a regime where Tj, is strongly sup- 
pressed with doping. In this regime, the s"* state is 
fragile and easily destroyed by the Coulomb repulsion, 
contrasting to the situation where all Fermi pockets are 
present. 

The paper is organized as follows: in Section II we 
discuss the competition between AFM and SC solely on 
phenomenological grounds. In Section III we introduce 
the mean-field Hamiltonian and derive the gap equations, 
the quasiparticle excitation spectrum and the free energy 
expansion. Section IV is devoted to the application of 
the formalism developed in Section III to the iron ar- 
senides, and is divided in 5 subsections. In Section IV-A, 
we present the results in the special case of a particle- 
hole symmetric band structure. Section IV-B contains 
both analytical and numerical results for various band 
structures without particle-hole symmetry. Phase dia- 
grams for parameters describing Ba (Fei_xCox)2 As2 are 
presented. In Section IV-C, we briefly discuss the regime 
where the sign of the coefficient that couples the AFM 



and SC order parameters becomes negative, and how it 
can be avoided by the onset of incommensurate AFM. 
Section IV-D discusses the effects of intraband interac- 
tions and Coulomb repulsion. In Section IV-E, we deter- 
mine analytically the shape of the reentrant Neel transi- 
tion line at low temperatures and the corrections due to 
fluctuations. In Section V we solve the same model pre- 
sented in Section III, but now with localized magnetic 
moments instead of itinerant AFM. Section VI discusses 
the suppression of T^ in the overdoped side of the pnic- 
tides phase diagram, and how it is related to the doping 
evolution of the Fermi surface. Section VII is devoted 
to the conclusions and, in Appendix A, we derive the 
Ginzburg-Landau coefficients in terms of Feynman dia- 
grams. Some of the results have been published in a short 
publication^^. 



II. PHENOMENOLOGICAL ANALYSIS 

Regardless of the microscopic details, the competi- 
tion between superconductivity and antiferromagnetism 
near their finite temperature phase transitions can be de- 
scribed in terms of a Ginzburg-Landau theory of coupled 
order parameters. The homogeneous part of the free en- 
ergy is given as 



Fgl (A, M) 



a 



V 2 



lAI 



^lAI^M^ 



-I— iiiM^ 



-M' 



(2) 



where A and M denote the SC and AFM order parame- 
ters, respectively. As usual, A is a complex order parame- 
ter, characterized by an amplitude and a phase, and M is 
a three component vector. The leading term in the order 
parameter competition is characterized by the coefficient 
7 > 0, where the sign of 7 reflects that both ordered 
states compete. As usual, the quadratic coefhcients are 
given by a™ = 0^,0 {T - Tn,q) and a^ = a^^o {T - T^fi) 
and change sign at T^.o and Tc.o, denoting the Neel 
and SC transition temperatures without order param- 
eter competition. We consider the situation where the 
transitions for 7 ~ are second order, i.e. the quartic 
coefficients m,„ and Us are positive. 

Furthermore, we consider that Tjv^o (2;) and Tc.o {x) 
vary as function of a physical parameter x that could 
be pressure, electron density or magnetic field. In case 
where both transitions meet at x = a::*, i.e. for 



rji^ 



Tn.o{x*) ^Tcfi{x*), 



(3) 



we have a multicritical point {x*,T*) in the phase di- 
agram. The vicinity of this multicritical point is the 
regime where a simultaneous expansion of the order pa- 
rameters is allowed. The mean-field analysis of Eql2] al- 
lows for two options for the phase diagram near {x*,T*), 
depending if 7^ > UmUg or 7^ < UmUg. Since we are 



interested in 7 > 0, it is convenient to define the dimen- 
sionless quantity^ 



7 



■y^llfJlllg 



- 1. 



(4) 



Thus, the nature of the phase diagram is determined 
solely by the quartic coefficients in the Ginzburg-Landau 
expansion Eql^J For g < Q (i.e. < 7 < y/UsUm), 
(x*,T*) is a tetracritical point where two second order 
phase lines cross, leading to a regime in the phase di- 
agram where simultaneous AFM and SC order occurs 
homogeneously within the sample, see Figl5^. In this 
regime both phases compete, but do not exclude each 
other. On the other hand, if 17 > (i.e 7 > y/u^u^) the 
phase competition is sufficiently strong that both phases 
are separated by a first order transition that terminates 
at the bicritical point {x*,T*). Notice that if the pa- 
rameter X jumps discontinuously from xi to X2 at the 
first order transition, there is an intermediate regime 
xi < a; < a;2 of heterogeneous phase coexistence, see 
FigHb- A sharp line of ffist order transitions occurs if 
one considers the phase diagram as function of h^, the 
variable that is thermodynamically conjugate to x, see 
Figl2t. Critical fluctuations, that go beyond this mean- 
field analysis, change the universal exponents near the 
critical temperatures and the slopes of the phase lines 
near (x*,r*). However, neither the generic behavior 
shown in Fig.l nor the quantitative criterion based on 
the sign of g are changed by fluctuations^!^. 

If we consider 5 < 0, both order parameters can be 
finite simultaneously. This regime is often referred to as 
coexistence of AFM and SC, referring to coexistence of 
order. This should not be confused with phase coexis- 
tence in the thermodynamic sense. The area in Figl^^ 
below the tetracritical point is a single thermodynamic 
phase characterized by two order parameters that are 
simultaneously finite. Similarly, the tetracritical point 
is not a point where four phases coexist (which would 
not be allowed by Gibbs phase rule), but a point where 
the system is in a single phase and both order parame- 
ters are infinitesimal simultaneously. Below the bicritical 
point, coexistence of thermodynamic phases only occurs 
for xi < X < X2, where macroscopic AFM and SC re- 
gions occur together in the sample. We use the term ho- 
mogeneous coexistence of AFM and SC order below the 
tetracritical point to refer to coexisting order and hetero- 
geneous coexistence below the bicritical point to refer to 
coexistence of phases. 

From the Ginzburg-Landau expression ([5]), we obtain 
the temperature dependence of the magnetic moment in 
the SC phase in the case of homogeneous coexistence: 



]y[2 /y^ ^ a-mfiUs (Tma - T) + 0^,07 (T - Tcfi) 
Y ^rn'^s 7 

Without phase competition both order parameters de- 
crease as function of temperature, dlsA^/dT < and 



T* 




T A 




(a) 



Xi Xj 

(b) 



sc 




Figure 2: Schematic phase diagrams {x, T) for competing AFM and SC orders. Here, a; is a generic physical parameter and T is 
the temperature. SoUd (dashed) lines denote second-order (first-order) phase transitions. For <? < 0, there is a tetracritical point 
and a region of homogeneous coexistence (a), whereas for g > Q there is a bicritical point (b and c). If a; changes discontinuously 
across the first order transition (panel b), forming a region of heterogeneous coexistence for x\ < x < X2 (shaded area), then 
its conjugate variable h^ changes continuously and the phase diagram has only one first-order line (panel c). 



dA/dT < 0. Phase competition can change this behav- 
ior. For instance, if 



as,Ql > arnfiUs, 



(6) 



it follows from Eq. (O that dM'^/dT > once su- 
perconductivity sets in. Thus, below Tc, the ordered 
moment decreases with decreasing temperature, as was 



observed in neutron diffraction experimenti 
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both 



Ba(Fei_xCox)2 As2 and^^ Ba(Fei_xRhx)2 As2. It is in- 
teresting that the same condition implies a back-bending 
of the antiferromagnetic phase boundary upon entering 
the superconducting state, see Fig. la. To demonstrate 
this we write for the bare Neel temperature Tn.o{x) = 
Tcfiil + f{x)), with df/dx < and f{x*) = 0. Thus, 
T/v,o is a monotonic decreasing function of x and meets 
the SC phase line at the carrier density or pressure 
value X* . Without restriction we assume that Tcfi is x- 
independent near x* . From the Ginzburg-Landau expan- 
sion Eq.(l2|) follows that the Neel temperature T/v inside 
the SC phase is given by: 



T, 



N 



Tc,o [l-fix]- 



amfiUs 



asfll - CLmfiUs 



(7) 



Since df/dx < 0, it follows that dT^/dx > in 
case Eq.(I5]) is valid. Thus close to x* one finds reen- 
trance of the paramagnetic phase below the SC transi- 
tion temperature. In Ba(Fei-xCox)2 As2 this was in- 
deed observediS^. The observation^^ of dyP/dT > in 
Ba (Fei_xRhx)2 As2 implies that the phase line for this 
material must bend back as well. Furthermore, the SC 
transition temperature inside the AFM phase is given by: 



Tc = Te,o 1 + / {x) 



am, 07 



(8) 



Using the condition Eq.(|n]) and the fact that 7 < 
y/UmUs {g < 0) it automatically holds that am,o7 < 



o-s,oUm, implying dTc/dx > 0. Hence, the SC transition 
temperature is suppressed inside the AFM phase. 

A very interesting limit is (7 = 0, i.e. at the transi- 
tion between the tetracritical point and bicritical point. 
Focusing on this multicritical point, where ag and «„ 
change sign simultaneously, one can introduce the five- 
component vector 



^=f^ 



1/4 



ReA, ImA,^M 



a„ 



(9) 



In case where the additional condition u^ = a^Um/cLm 
is fulfilled it follows that the free-energy can be written 
as: 



FGLiA,M)^Jd''r[^f^' + ^^' 



(10) 



where a — sign(as) \/\as\ \am\ and u = ^JUsUm- This is 
the SO (5) symmetric form of the Ginzburg-Landau en- 
ergy that was first proposed by Zhang to describe the 
physics of the cuprate superconductor a^'*'^^ . In the con- 
text of the pnictides, it was shown^ that a model Hamil- 
tonian similar to the one used in this paper is invariant 
with respect to a global SO (6) symmetry that contains, in 
addition to the AFM and SC order parameters, an imag- 
inary density wave state. Below we will see that there 
is evidence that the the pnictides are indeed strongly af- 
fected by such an enhanced symmetry. 



III. MICROSCOPIC MODEL 

So far we have analyzed the problem of coexistence be- 
tween SC and AFM only on phenomenological grounds, 
which gave us interesting and general information about 
the phase diagram. Next we develop a microscopic model 
that captures the essential aspects of the iron arsenides 
to determine their detailed phase diagram. We will also 



make explicit contact to the Ginzburg-Landau theory 
and determine the coefficients of the order parameter ex- 
pansion to obtain the behavior close to the transition 
temperatures. 

We start from the Haniiltonian: 



v. —Ho + Hafm + "Hsc- 



(11) 



The non-interacting part Ho describes two bands 
shifted by the momentum Q relative to each other: 



Ho 



ko 



(ei,k - m) cLck^+E ('^2,k+Q - fJ.) 4+Qo-'^k+Q^- 



ko- 



(12) 

We consider only one hole band located in the center of 
the Brillouin zone with dispersion ei,k, and one electron 
band, shifted by Q from the hole band, with dispersion 
£2.k- To keep the discussion as simple as possible, but 
still capturing the basic properties of these materials, we 
consider a circular hole band and an elliptical electron 
band: 



ei,k 

-2,k+Q 



— £1.0 ^ 



2™' 



-£2,0 + 



2m^ 



2m,, 



(13) 



where ea,o is the energy offset (see figure [J]). Such 
a choice is motivated by angle-resolved photoemission 
spectroscopy™ (ARPES) as well as by tight-binding fit- 
tings to first-principle band structure calculations. The 
operator cj^^ ('^k+Otr) creates an electron with momen- 
tum k (k -H Q) and spin a in the hole (electron) band. 
The chemical potential is denoted by /i, and we define 
^Q,k = £Q.k~/^- Frequently, it will be possible to gain ex- 
plicit analytic insight by considering the limit of particle- 
hole symmetry, where Eq = £1,0 = £2,o> 'm^x = iTT-y = m 
and 11 — 0. In this case the hole and electron Fermi 
surfaces are identical, leading to perfect nesting. In what 
follows we supplement our numerical analysis of the prob- 
lem with band structure Eq. p^ by analytical results at 
or near the limit of particle-hole symmetry. Even though 
all five iron d-orbitals contribute to the states at the 
Fermi surface of the iron arsenides, the physics of the 
competition between the antiferromagnetic and the su- 
perconducting states is well captured by this effective 
two-band modeP°. 

The same electrons that form the Fermi surface are as- 
sumed to be responsible for the magnetism of the system 
through an electronic interband Coulomb interaction J, 
leading to excitonic itinerant antiferromagnetism^SiSli^: 



HaFM = I 2J E'^ks'^ss'rfk+qs' •dk's'^ss'Ck'-qs'- (14) 
k,k',q s,s' 

Here ct^^, denotes the {ss') element of the i-th Pauli 
matrix, with s = ±1. In the weak-coupling limit, we 




Figure 3: Schematic representation of the two-band structure 
used here. It consists of a circular hole band (band 1, red line) 
in the center of the Brillouin zone and an elliptical electron 
band (band 2, blue line) shifted by the ordering vector Q 
from the hole band. The chemical potential /i is positive for 
electron doping. 



perform a Hartree-Fock decoupling that leads, except for 
a constant shift in energy, to the effective single particle 
Haniiltonian 



Ha 



FM 



ks 



sM 



^ks^k+Qs 



^k+Qs 



Cki 



(15) 



where N is the system size. M denotes the antiferromag- 
netic gap opened at momenta ko that are Bragg scattered 
due to magnetic order, ^i,ko = '^2,ko+Q- In general, for 
small M, the Fermi surface is only partially gapped and 
the magnetic state is metallic. For large enough M, how- 
ever, the Fermi surface can become completely gapped. 
In the case of perfect nested bands, an infinitesimal anti- 
ferromagnetic gap is already able to gap the entire Fermi 
surface. Note that M — |M| is proportional to the am- 
plitude of the staggered magnetization m and given by: 



k,s 



/ m. 



(16) 



Besides the magnetic interaction /. the electrons arc 
also subject to a pairing interaction Vap, where a, P — 
1, 2 are band indices. In case of pure interband interac- 
tion Vap = V {1 — Sap) the Hamiltonian becomes 



Hsc = V 



-k+qt''-ki^ 



^tC^k't- 



(17) 



k,k' 



Below we demonstrate that the introduction of an in- 
traband pairing interaction does not change the main 
conclusions of our paper. A„ is the superconducting 
gap of band a which, given the interband coupling V, 
is due to the action of the electrons in the opposite band 
a. Thus, Ai and A2 are determined by the two gap- 
equations: 



V 



^^ - -i/E V- 



N 
V 



-ki 



^1 = -i^E 4+Qt^V 



Qi 



(18) 



The expectation values are determined with the mean 
field Hamiltonian: 



H 



sc 



E(^i4/-H + H (19) 

k 

E (^24+Qt'^-k-Qi + '^-C- 



The mean-field Hamiltonian formed by the sum of Eqs. 
(fT2| . p3|) and (fT9| is quadratic and can be diagonalized 
analytically, yielding the self-consistent gap equations: 



An 



^^i^A^tanh^^^- 



TV 



M = ly^ivrf^tanhf^^ 

^ Z^ k,a I 2 

1, ^ \ 



(20) 



with kernels: 



K^ - 






A, 



a (-^aJ 



A2 _ f2 



M^An 



2Ea 



k - ^' k 
,k a,k 



M 



(£;2k + AiA2 + Ci,k6,k-Af' 



2^;, 



'a^k (i?i 



k -^a.k 



(21) 



The excitation energies Ea^k (a = 1, 2) are the positive 
eigenvalues of a state with simultaneous magnetic and 
superconducting order: 



E' 



rw±A/r: 



rjL 



(22) 



(1^0)1 follow as the stationary points df (M, Aq.) /dAa — 
df{M,A^)/dM = 0. 

Before we analyze the impact of magnetic long range 
order on the pairing state we discuss the gap equations 
(|20p in the limit M = 0. Here, we perform the momen- 
tum integration by introducing the density of states pi 
and p2 of the two bands. The condition for T^ is that the 
largest eigenvalue of 



A = 



-Vp2 

-Vpi 



(24) 



is positive and equals to 1/ \n{W/ (aTc)), were a 



'lE 



/2 , ^E is Euler's constant and W is an upper 
energy cutoff for the pairing interaction. Clearly, the 
eigenvalues of A are A± = ±.V^/pip2- The pairing state 
is determined by the corresponding eigenvector 



(Ai,A2)cx 



1 



VP1+P2 



(\/P2, T^/Pl) 



(25) 



Thus, for F < 0, A_ is the largest eigenvalue and the 
gap equation has a solution where Ai and A2 on the 
two sheets of the Fermi surface have the same sign. This 
is called the s++ state and is analogous to the pairing 
state of the multiband superconductor^^ MgB2. It is 
the natural state that arises as a result of conventional 
electron-phonon interactions. On the other hand, the 
gap equations also allow for a solution for V^ > 0, with 
Ai and A2 having different signs on distinct Fermi sur- 
face sheets. This s"' -state results from purely electronic 
interactionsr^i^"— 1^"— i.e. it is the natural analog to 
the (ia;2_j,2 -pairing state in the cuprates with a single 
Fermi surface sheet. As the relative sign of Ai and A2 
is —1, this is an unconventional SC state, even though 
it is in the same irreducible representation Aig of the 
symmetry group D^h as the s"'""'"-state. 

In case where Ai and A2 have the same sign (s++- 
state), the excitation energies ±i?2,k possibly have 
nodes-^*'. The nodes are located at the set of points k„ 
that satisfy simultaneously the equations: 



A2 
^2 



'k,l 



+ e- 



k+Q,2 



and 



with Fk = 2M-^ + Af 

rik = -4 fe^,! + ^1) (4+Q.2 + ^2) as well as (5k = 

8M2(AiA2-f £k,iek+Q,2-MV2). The free energy 
density of a system with SC and AFM long range order 
that results from this analysis is: 

/(M,A„) = ^M2-l(AtA2 + A;Ai) 



— >^ log I 2 cosh I — — — 



.(23) 



The superconducting order parameters Ai and A2 
of the two bands and the staggered moment oc M are 
obtained by minimizing /(M, A^). The gap equations 



;i,k„ 



6.k„+Q = ±./|^(A/2-AiA2)'/' (26) 

V ^1 

determined from the condition i?2,k = 0. Obviously, one 
condition for nodes to exist is Af^ > A1A2. For Ai — A2 
the condition for nodes in the antiferromagnetic state cor- 
responds to ^i.k„ = C2.k„+Q i-e. where Bragg scattering 
due to antiferromagnetism is large. However, nodes are 
not guaranteed to emerge. For example, in the case of 
particle hole symmetry (implying perfect nesting of the 
Fermi surface) 



a^e 



l.k 



-e 



2,k+Q 



(27) 



it holds for / > |V| that A = Ai = iAj = -l^M < M 
and the above equations cannot be fulfilled simultane- 
ously. This also follows if one explicitly considers the 
eigenvalues for s++ pairing in the limit of particle-hole 
symmetry: 



■E^a.k 



^l 



(M ± A) 



(28) 



which are fully gapped for Af ^ A. Only for |y| = / fol- 
lows an entirely gapless eigenvalue i?2,k, consistent with 
the condition Eq. ipS)) . In distinction, the eigenvalues for 
s"' pairing and particle-hole symmetry are fully gapped 
and doubly degenerate: 



■£^a,k - ^k 



Ar 



(29) 



It is also interesting to note that, if one considers the 
simplification A = Ai = — A2 even in the absence of 
particle-hole symmetry, then the excitation energies for 

the s"' case assume the simple form -Ef k ~ ( -^a\ 
A^, where 



.r^AFMV 



E. 



AFM _ /Cl,k + 6.k 



a.k 



±\ AP 



^l,k — ^2,k 



(30) 



are the excitation energies of the pure AFM state. Thus, 
in this special situation, one can perform two sepa- 
rate Bogolyubov transformations to diagonalize the full 
Haniiltonian. 

These considerations allow for some general conclu- 
sions of the order parameter dependence of the free en- 
ergy Ea. ((23|) . The quadratic terms of / (Af, A^) depend 
on the interaction strengths / and \V\. However, all other 
dependencies take place only via the implicit dependence 
of Ea,k on the order parameters. In case of s++ pair- 
ing, the order parameters enter these extra dependencies 
through the combinations (Af ± A) , while for s^ pair- 
ing the third term in Eq. (P5|) can only depend on the 
combination Af^ -I- A^. Thus, we find for the free energy 
of the s^ -state with particle-hole symmetry: 



/+- = 



2Af^ 



2A2 
M 



2$ Af 



(31) 



On the other hand, it follows for the s++-state: 

- -^^^^2^^<I>f(lA/|+|A|)^ 



/++(A,Af2) = _ 



11^1 



$ (|Af|-|A|)^ , 



where 



/•OO 

$ (x) = -4Tp / dz 

J\x\ 



zlog(2cosh(/3z/2)) 
\/z^ — x^ 



(32) 



(33) 



is the same function in both cases. These facts have 
important implications for the Landau expansion of the 
free energy that we discuss next. 



In order to obtain microscopic expressions for the 
coefficients of the Ginzburg-Landau theory, we expand 
(5/(M,A„) = /(M,A„) - /(0,0) with respect to M 
and Aq . It follows from Eq. (|23p that 

6f = ^M^ + ^M^ + ^^A^A^ 

+ 5:^Al+5:2|^M^A„A, (34) 

For the coefficients of the antiferromagiietic order pa- 
rameter follows: 



am = J - 2xph (Q) (35) 

1 ^i.ksech^ (^) - A2,k+Qsech2 (%tQ) 



N ^^ 

k 



7"(Cl,k - C2,k+Q) 



with coefficients ^Q,k = —£,a,k + £,a,u+Q + 
2Tsinh(^Q_k/T) and the bare static particle-hole 
response at the antiferromagnetic ordering vector: 



Xph 



(Q) = ^E 



l^tanh(%^)-tanh(%Q) 



N . Cl,k — 'C2,k+Q 



(36) 



For the coefficients of the superconducting order pa- 
rameters follows : 



as,a/3 = -y {i- Sap) - SapXpp{0) 



(37) 



ANT 



E 



sech^ (%ii) (Tsinh {/3Um) " C«,k) 



'5Q,k 



where 



-^ tanh 

k.a 



2T J 



(38) 



is the bare static particle-particle response at external 
momentum q = 0. Finally, it follows for the coefficients 
7q/3 that determine the coupling between both order pa- 
rameters: 



1 ^ 



c. 



ii) 



^"" 27V 2^^ TC^^ (^„^k + ^, ,) (^„^k _ ^,^^f ^^^) 

_ J_ v^ Cojc tanh (/3Ca,k/2) - ^a,k tanh (^^a,k/2) 



.(1) _ 



J^ ?a,kSa,k ( Sa^k ^a.\ 



with C^^i =. 2T tanh (/3^a^k/2) ^^,k - C^k + ^UkUk , 



C 



S - -ea.k(elk-€,k)sech^(/3ea,k/2) as well as 

.(3) 



Q:k--4T^o,ktanh(/3^a,k/2). 
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Figure 4: Diagrammatic representation of the quartic 
Ginzburg-Landau coefficients associated to tlie system with 
competing AFM and SC order parameters. The single- 
particle Green's functions of the two bands are denoted by 
Gi (k). 



These Ginzburg-Landau coefficients can also be ex- 
pressed in terms of Feynman diagrams, obtained by in- 
tegrating out the fermionic degrees of freedom of the 
system with competing AFM and SC. The derivation 
is presented in Appendix A; in Fig. |4l we show the 
diagrammatic representation of all the quartic coeffi- 
cients in terms of the single-particle Green's functions 
Gi (k) = (jw„ -(zmV ■ 

Due to the coupling between the two bands, Ai and 
A2 will always appear simultaneously. As follows from 
the eigenvectors of A in Ea. ([M| . close to Tc, the ratio 

A1/A2 — ±(p2/pi) is determined by the ratio be- 
tween the densities of states of the two bands. In our 
case holds P2I P\ = JrHxiTiy/m. The relative sign of Ai 
and A2 depends on the sign of V . Thus, one can intro- 
duce the superconducting order parameter A via 



Ar = J^A 

Pi +P2 



A2 = ±W^A 

Pi +P2 



(40) 



such that A^ = (Af + A|) /2. This leads to the Lan- 
dau expansion coefficients of the superconducting order 
parameter of Eql2] 

2as.iiP2 + 2as^22Pi + 4 |as,i2| VP1P2 

"« = ' 1 \ ^ 

(P1+P2) 

4us,ip| -I- 4ms 2pi 

Us — 2 

(Pl + P2) 

2711P2 + 2722P1 ± 47121/pIpJ 
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(41) 



(Pi + P2) 
Note that the coefficient 7 of the coupling between the 



SC and AFM order parameters depends on the relative 
phase of the two SC order parameters. In the next section 
we analyze these expressions and discuss the implications 
of these results for the phase diagram of the pnictides. 

Before moving on, let us briefly discuss the relation- 
ship between the system's dimensionality and the Neel 
transition temperature. Since our mean-field model is 
insensitive to fluctuations, it allows a flnite T^ even for 
two dimensions. However, due to Hohenberg-Mermin- 
Wagner theoreniS, a non-zero Tm will only be possible 
if the band structure has a three-dimensional dispersion, 
i.e. if the electronic interaction responsible for the AFM 
instability is effectively 31?. In the iron arsenides, the 
AFM ordering involves Fe ions located on spatially sep- 
arated layers. Assuming a weak interlayer coupling, we 
can introduce a phenomenological 31? anisotropic action 
for the low-energy collective magnetic modes and obtain 
the Neel transition temperature T^r: 



T, 



N 



Tn 



^^-•mr^^ 



(42) 



Here, J is the effective in-plane magnetic exchange, J^ 
is the interlayer coupling and T^ is the mean-field Neel 
transition temperature. Notice that J and J^ are effec- 
tive parameters of the low-energy model originated by the 
electronic interaction Eq. p^ . and are not necessarily 
related to localized spins. The logarithmic dependence of 
T/v on Jz/ J shows that the overall scale of the transition 
temperature of an anisotropic magnetic material is deter- 
mined by the mean-field value Tj^. This explains why, in 
the iron arsenides, Tm has the same order of magnitude 
for both the 1111 and the 122 compounds, even though 
the former are much more anisotropic than the latter. 



IV. PHASE DIAGRAMS 

In this section, we will use the formalism developed 
above to explore the possible phase diagrams of the sys- 
tem with competing SC and AFM order. In particular, 
we will be interested in analyzing whether different su- 
perconducting states are able to coexist with an itiner- 
ant antiferromagnetic state or destined to phase separate 
from it. 



A. Particle-hole symmetric case 

As we stated in Section III, the band structure of the 
iron arsenides can be generically described by two sets of 
hole and electron bands, displaced from each other by the 
magnetic ordering vector Q. Even though the electron 
and hole bands are not perfectly symmetric to each other, 
we can start our analysis by considering, at first, the case 



of two nested bands, such that ^k = Ci k 



-6, 



k+Q- 



No- 



tice that, in this context, nesting does not mean that the 



distinct pockets of the Fermi surface have parahel seg- 
ments; instead, it imphes that they have the same shape 
and area, such that the non-interacting Hamihonian Hq 
has particle-hole symmetry. 

The limit of perfect nesting corresponds to Eq = Si^q — 
£2,0 J 'Tij; = my — m and fi ^ in Ea. (jl3p . In this case, 
it is straightforward to conclude that A — |Ai| — IA2I. 
Moreover, using formulas ([37]), it follows that Og^n = 
as,22 and Us_i — Us,2- Thus, independently of the rela- 
tive sign between the Cooper pair wave functions of the 
two bands, they have the same SC Ginzburg-Landau co- 
efficients, meaning that the thermodynamic properties of 
the "pure" s~^^ and s^ states will be the same. How- 
ever, the coupling to the magnetic degrees of freedom 
significantly changes this picture. 

The Ginzburg-Landau expansion is formally valid 
around the temperature where the AFM and SC phase 
lines meet, Tn — T^. From the magnetic and supercon- 
ducting quadratic coefficients, we see that, for particle- 
hole symmetry, this condition implies I — \V\. Thus, cal- 
culation of the Ginzburg-Landau coefficients using Eqs. 
dSSl), dSZl) and dMl) yields: 



4 2 



(X — Clm, ' ^,s T 



4: Z. ^— V 



tanh 



\^t) 



1 (W 



where a = Tre ^'^ /2, as well as: 



(43) 



Urn = Us = 711 = 722 = 27; 
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(44) 



with: 



2NT 



E' 



scch" 



(|f.)(rsinh(^) 



-6 



^3 



P 7C(3) 



(45) 



Inserting these results into the Ginzburg-Landau ex- 
pansion Eq. (p4|) yields^S: 



<5/(M,A) = -(M^ + A^) 



4 ^ 



A^ 



2 a2 



^ 2 



(46) 



where g = {1 + cos6'), as given by Eq. (j4|), with 6 denot- 
ing the relative phase between the SC order parameters 
of the electron and hole bands. Thus, for the s++ state 
{0 = 0), it follows that g = 2 > 0, meaning that the 
s"*""*" state is deep in the mutual exclusion regime, unable 
to coexist with AFM in the region of the phase diagram 
close to Ttv — Tc- However, for the s"' state it holds 



that 5 = 0, implying that this state is in the borderline 
between the coexistence and mutual exclusion regimes. 

We can trace back to the Ginzburg-Landau coefficients 
in Eqs. ([55]) . ([57)1 and (15^ the origin for the distinct be- 
haviors of the systems with competing itinerant antifer- 
romagnetism and s^^ or s^ superconductivity. As we 
showed above, the quadratic and quartic SC and AFM 
coefficients are the same in both cases. However, the net 

SC-AFM coupling coefficient 7^ = 711 -I- 722 — 2712 is 

reduced in the case of an s"* state when compared to 
the case of an s~^~^ state, where 7++ = 711 + 722 + 2712. 
Notice that in both situations 7 > 0, evidencing the com- 
petition between the two phases. 

In fact, from the diagrammatic representation of the 
coefficients presented in Fig. (H it is clear that the 
only Feynman diagram sensitive to the relative phase 
between the SC order parameters of the two bands is 
the diagram corresponding to 712 . It gives a contri- 
bution to the free energy of the form M^ |^i| |^2| cos6', 
see EqUJ Therefore, the static long-range magnetic or- 
der plays the role of an intrinsic Josephson coupling: it 
provides the momentum Q to the electrons of a Cooper 
pair in band 1, scattering them coherently to band 2, 
where they recombine. Thus, the region of the phase di- 
agram where antiferromagnetism and superconductivity 
compete provides an efficient tool to probe the relative 
phase between the Cooper pair wave functions, an in- 
formation that is usually reserved to very delicate phase 
sensitive experiments^"—. The existence of such a tool 
is even more relevant in the case of the iron arsenides, 
since both s++ and s^ states belong to the same irre- 
ducible representation Aig of the tetragonal point group 
Dih, making interference experiments rather involved 
and comple x^^i^^" — . 

The analysis of this limiting case with particle-hole 
symmetry suggests that while the s"'""'" state is intrinsi- 
cally unsuitable for coexistence with the AFM phase in 
the iron arsenides, the s"' state may or may not coexist 
with magnetism. In the realistic case where particle-hole 
symmetry does not hold, the decision on whether the 
s^ state is in the regime of coexistence or mutual 
exclusion will depend on additional details of the band 
structure, as we will demonstrate in the next subsection. 
This explains why some compounds present homoge- 
neous coexistence^i^iii^"— , like Ba(Fei_xCox)2 As2 
and possibly^ SmFeAs (Oi_xFx), while in oth- 
ers, such as LaFeAs(Oi_xFx), PrFeAs (Oi_xFx), 
(Sri„xNax)Fe2As2, and (Bai_xKx) Fe2As2, AFM and 
SC are mutually excluding^i^rJ^. 

Notice that these results do not depend on the spe- 
cific functional form of the bands dispersion relations 
nor on the dimensionality of the system. They follow 
solely from the assumption of particle-hole symmetry 
^i.k = — ^2,k-i-Q- Note also that, as we anticipated in 
the previous sections, the free energy ([^5[) for the s^ 
case is completely symmetric with respect to both or- 
der parameters and can be characterized by the SO (5) 
order parameter N = (Re A, ImA, M). Remarkably, a 
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similar S0(5) model has been proposed previously for 
the cuprates5M5. In the context of the iron arsenides, 
recent works^^i^i demonstrated that the complete, inter- 
acting particle-hole symmetric Hamiltonian has an emer- 
gent SO (6) symmetry (the other degree of freedom which 
is not captured in our model is associated to an imagi- 
nary density wave). The existence of such an emergent 
symmetry suggests that our result regarding the ability 
of the s^ state to coexist with magnetism is likely valid 
not only in our weak-coupling approach but also in the 
strong-coupling limit. 



B. General case and application to 

Ba (Fei_xCox)2 As2 

We now move away from the particle-hole symmet- 
ric case and consider more specific details of the band 
structure dispersions of the iron arsenides. Let us first 
consider small perturbations that break the particle-hole 
symmetry. For instance, we first take Eq = £1.0 = £2.0j 
nix = iTT-y = m, but a finite chemical potential ^ <^ T, 
i.e. we have two detuned circular bands. An analytic 
expansion yields, to leading order: 



fl (meV) 



g+- w 0.018 
.9++ 



M 



T 
2 + 0.386 f^ 



(47) 



Thus, in the case of spherical detuned bands, we al- 
ways find a first-order transition between the supercon- 
ducting and magnetic phases, independent of the pairing 
state. This is in agreement with numerical calculations 
performed by Vorontsov et ai, which found no region of 
coexistence between commensurate AFM and SC"^*. 

The second perturbation we consider is an infinitesimal 
ellipticity of the electron band, such that rrix — m + Sm 
and my — m — dm, but with eg = £1,0 — £2,0 a-nd /i — 
0. Such perturbation also makes |Ai| and IA2I assume 
different values. In this case, we obtain the following 
perturbative expansion of g for the s"' case: 



9+- 



-0.0039 



0.0022 ( — 



-F0.00008 






6m 
m 



(48) 



while (7-|_+ remains close to 2. Since £0 3> T, we conclude 
that the s^ state moves again to the regime of mutual 
exclusion from antiferromagnetism. 

In order for the s^ state to be able to coexist with 
AFM, we need to consider both a finite chemical potential 
and a finite ellipticity. Then, depending on the particu- 
lar values of the band masses, of the energy offsets and 

of the chemical potential, g^^ will be either positive or 

negative, while (7++ remains positive. 
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Figure 5: (a) Ellipticity of the electron band e = , /l — ^ 

and (b) chemical potential \jl (in meV) as functions of t. The 
parameter t interpolates between two points of the band struc- 
ture parameters space: i = corresponds to the particle-hole 
symmetric case while t = \ refers to the parameters that give 
good agreement with experimental magnetization data^" on 
Ba (Fei_xCox)2 As2 (see the main text for the actual values). 



To illustrate this, we perform a numerical calcula- 
tion of the coefficients g through a particular path con- 
necting two important points of the parameters space. 
They are the point with particle-hole symmetry, where 
£0 = £1.0 — £2,0i "^x = ^^y = "^ a-nd /i = 0, and the point 
corresponding to the values that consistently describe the 
magnetic properties of Ba(Fei_xCox)2 AS2: £0 = 0.110 
eV, £1,0 = £0 - A^o, £2,0 = £0 + Mo, Mo = 0.015 eV, m = 



1.32m( 



Icctr 



Im, TO,, 



0.3?7i and 11 — fie = 0.039 



eV. As we showed in a previous work^S, these parameters 
give a satisfactory agreement between our model and the 
doping and Tjv dependence of the experimental values of 
the relative zero-temperature staggered magnetization in 
the absence of SC, M {x,T ^ 0) /M {x = 0,T ^ 0). The 
chemical potential /Xc corresponds to a variation of the 
electronic occupation number by An w 0.06 with respect 
to /iQ. Since each added Co atom replaces one Fe atom, 
we associate this increase of An to the Co doping con- 
centration X = 0.06. 

In particular, the path chosen to connect these two 
points is parametrized by a real number t G [0,1], 
such that TOj; = (l-l-i^jTO, niy = {1 — 0.7 t)m and 
/i = (/ic + /io)i^- The variation of the electron band 
ellipticity and of the chemical potential as function of 
t is shown in figure [SJ and the corresponding values of 

g^ and <;++ are presented in figure [6] Clearly, when 

both ellipticity and fj, are finite, g^ can be either posi- 
tive or negative, but g+^ remains positive. Notice that, 
for the parameters corresponding to Ba (Fei_xCox)2 As2 
{t = 1), the s^ state coexists with magnetism, while 
the S++ state is incompatible with AFM. In particular, 
for t = 1, we have: 



9+- 

9++ 



-0.52 
2.0 



(49) 



Our analytical results and numerical calculations indi- 
cate that (7-(_+ is generally positive, in special for the range 
of parameters associated to the pnictide compounds. 
In order to investigate this point further, we analyzed 
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Figure 6: Coupling coefficient g for both the s"' state (a) and 
the 3+"*" state (b) as function of the parameter t, which inter- 
polates between distinct band structure parameters, changing 
simultaneously the electron band ellipticity and the chemical 
potential (see Fig. [5] for the definition of t). 



which parameters are able to bring g_|_+ to smaller val- 
ues. Particularly, we noticed that by increasing the mass 
anisotropy of the electron band, while keeping the chem- 
ical potential fixed, g++ can be reduced. In figure [7l we 
show the effects of an extremely large mass anisotropy on 
the value of 5++. All the other band structure parame- 
ters have the values used before for i = 1. Clearly, even 
after pushing the electron band ellipticity to unphysical 
limits - at least in what concerns the iron arsenides - we 
still obtain that the §+"'' state cannot coexist with mag- 
netism. 

These results were briefly discussed in our previous 
work^° and, in some detail, by Vavilov et alM: Con- 
sidering that £1,0, £2,0 s-re the dominant energy scales of 
the problem, the authors of Ref.'*^ write the band struc- 
ture ()13p in the form ^2,k+Q = — ^i,k — 2(5^,, where ip is the 
angle along the electron pocket and S,p = Sq + 62 cos 2ip, 
with 60 proportional to the chemical potential and band 
masses difference and 62 to the ellipticity. In this limit, 
one can approximate |Ai| « IA2I and expand in powers 

of 60 and ^2. They obtain that g.^ becomes negative 

for a significant range of values where both ^o and 62 are 
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Figure 7: Coupling coefficient g for the s+"^ state as function 
of the electron band anisotropy mx/my. The values of the 
other band structure parameters are described in the text. 



simultaneously finite. Furthermore, they also find that 
(7++ is always positive, in complete agreement with our 
previous and present results. 

Although the Ginzburg-Landau expansion is extremely 
useful to investigate if SC and AFM are able to coex- 
ist, it is formally not valid far from the point where 
the two phase lines meet. In order to obtain a com- 
plete (a;, T) phase diagram, including the back-bending 
of the Neel transition line predicted phenomenologically 
in Section II, we self-consistently solve the gap equations 
pop at a fixed occupation number. Using the parame- 
ters discussed above for Ba (Fei_xCox)2 As2, we obtain 
the phase diagrams presented in figure [TJ A zoom of the 
phase diagram associated to the s^ SC state is pre- 
sented in figure [8j evidencing the reentrance of the AFM 
transition line. The magnitudes of the electronic inter- 
actions were chosen to yield^S T/v — 140 K at a; = and 
Tc = 25 K at a; = 0.06, and are given by \V\ = 0.46 
eV and / = 0.95 eV. The level of Co doping x is associ- 
ated to the variation of the electronic occupation number, 
which depends on the chemical potential. Specifically, we 
consider that each added Co corresponds to one electron 
added in the system. 

We emphasize that all band structure parameters were 
determined in our previous work^'' by fitting the Tn and 
X dependence of the experimental zero-temperature mag- 
netization, M {x, T = 0) /M{x = 0,T = 0), in the ab- 
sence of SC. Therefore, in the phase diagram presented 
in figure [U all the available free parameters are fixed by 
the shape of the transition lines Tn^q (x) and Tcfi{x) of 
the independent, uncoupled phases. The actual transi- 
tion lines Tn (x) and Tc (x) of the system with coupled 
AFM-SC phases are the solution of the self-consistent 
gap equations, with no extra free parameters involved. 

Clearly, the only difference between the phase dia- 
grams for an s"'""'" and an s"' SC state is on the coex- 
istence / mutual exclusion regions. In figure |9l we com- 
pare the temperature dependence of the AFM and SC 
gaps for a fixed doping in both cases. For T < T^ while 



12 



50 




S state 



SC 



(a) 



0.04 0.05 0.06 0.07 0.08 



Figure 8: Zoom of the phase diagram of Ba (Fei-xCox)2 As2 
considering a superconducting s"' state (see figure [l]). Note 
the reentrance of the non-magnetically ordered phase at low 
temperatures. 



in the s"' case the magnetization is strongly suppressed 
but still survives, in the s"^"*" case it completely vanishes 
once the SC gap opens. Note that, for the s"' case, the 
T = values for M and Aq are smaller than their values 
in the respective pure states. In figure [TUl we present the 
temperature dependence of the magnetization for several 
doping values in the case of AFM competing with s"* 
SC, demonstrating its stronger suppression as the tetra- 
critical point is approached. 

Numerical calculations of the phase diagram associated 
to the simplified band structure <^2,k+Q = — Ci,k ~ 2(5,^ 
discussed above were recently presented by Vorontsov et 
al^. Our results from figure [1] are in general agreement 
with their findings. Exploring other regions of the pa- 
rameters space, they also found systems where the s^ 
coexistence region does not persist all the way to T = 
as well as a small region at very low temperatures where 
s'^'^ could in principle coexist with AFM. 

A rather small region with coexistence between 
isotropic s-wave SC and itinerant AFM was also found by 
Kato and Machida'^'* in the context of heavy fermion com- 
pounds (see Section V for a brief discussion about these 
materials). Considering a single band without particle- 
hole symmetry, they performed numerical calculations to 
determine the phase diagram for different pairing states. 
In particular, coexistence between isotropic s-wave SC 
and AFM was only found far from the multicritical point 
T/v — Tc and in a very narrow regime, analogous to what 
was reported by Vorontsov et alA^ in the context of the 
iron arsenides. Note that these results are not in contra- 
diction to our conclusions, since our Ginzburg-Landau 
expansion - and, consequently, the definition of the cou- 
pling parameter g^^ - is only valid for T^ — Tc- Far from 
the multicritical point and from particle-hole symmetry, 
the details of the bands dispersions are very important 
and it is in principle possible to find coexistence even if 
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Figure 9: Absolute values of the superconducting order pa- 
rameters Ai and A2 (in meV), as well as of the magnetic order 
parameter M ((in units of its value Mo at {x = 0,T — 0))) 
as function of temperature T (in K) for the fixed doping 
level X = 0.054 (see the {x,T) phase diagrams of figure [l]). 
Panel (a) shows the result for an s"' state, whereas panel 
(b) presents the result corresponding to an s^^ state. 



.g++ > at Tat ~ Tc. 

So far, we have only compared the s^ and s++ SC 
states in our calculations. However, electronic theories 
for the superconductivity in the iron arsenides have also 
proposed other symmetries for the Cooper pair wave 
function where nodes are present^2iZ^i2&. One example is 
the d-wave state, where Ai and/or A2 have nodes along 
their respective Fermi pockets. The generalization of our 
formalism to these other symmetry states is straightfor- 
ward. One has only to introduce the corresponding an- 
gular factors 77 (if) for the gaps and for the pairing inter- 
action V , and then average over the Fermi pockets. 

In the case of particle-hole symmetry, where the free 
energy is given by Eq. (I46[) . we obtain, for a d-wave 

state, g = (\/f~l)'"~ 0.6u. Unlike the s^ state, the 
d-wave state is not on the borderline between coexistence 
and mutual exclusion from AFM. However, it is neither 
deep in the mutual exclusion regime, as the s'^'^ state 
is. Even though the s"' state is the most compatible 
with itinerant magnetism, we cannot exclude that the d- 
wave state is also able to coexist with AFM for certain 
parameters. 
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Figure 10: Magnetic order parameter M (in units of its value 
Mo at (x = 0, T = 0)) as function of temperature T (in units 
of the AFM transition temperature TV) for different doping 
levels. Due to the competition and coexistence with the s^^ 
SC state, M decreases below T^. Note the reentrant behavior 
for X = 0.059. The right panel is a zoom of some of the curves 
from the left panel. 



C. Incommensurability and the sign of the 
coupling coefficient 

Our model assumes that the magnetism is commen- 
surate. Experimentally, this is still an unsettled issue 
for Ba (Fei_xCox)2 As2: while neutron diffraction mea- 
surements did not detect any incommensurability inside 
their resolution windo w^^'^^i^'' , some works employing 
NMR^'^ and Mossbauer spectroscopj*21 suggest that the 
magnetism could be weakly incommensurate in these sys- 
tems. 

Theoretically, the weak-coupling model for the ex- 
citonic itinerant magnetism naturally predicts the on- 
set of an incommensurate AFM state for small enough 
temperatures^, as recently pointed out by Vorontsov et 
a0^. To see how this comes out from the model we used 
here, consider the specific case of detuned bands having 
the same shape, ^2,k+Q = ~Ci.k ~ 2/i. Instead of ex- 
pressing the Ginzburg-Landau coefficients as momentum 
sums, Eqs. pSI) . ([57| and (|39p . we can equivalently ex- 
press them as Matsubara sums. Using the diagrammatic 
form of the coefficients (see figure |3] and Appendix A) , 
this is a straightforward calculation and yields: 



47:pTj2 



UJn {l^l - 3^^) 



^„>o K + M^)^ 



2^pT y 4 



<^„>oKHm^)^ 



w„>0 



UJn {ujI + p?) 



(50) 



Similar expressions were obtained in Ref4i. Using Eq. 
([SU]) . it becomes clear now that Um < for T^ < Q.bfi, 
indicating that the transition from the paramagnetic to 
the commensurate AFM phase is first order. However, 
as discussed elsewhere^, when this condition is met an 
incommensurate AFM phase has a lower energy than the 
commensurate state. For the parameters we used to ob- 
tain the phase diagrams of Ba (Fei_xCox)2 As2 (Fig. [1]), 
where the Fermi pockets have actually different shapes, 
the AFM phase line meets the SC phase line before this 
incommensurate instability takes place. Even if they met 
after this instability point, it was shown by Vorontsov et 
al. that the s"* state remains able to coexist with an 
incommensurate antiferromagnetic state^Sii^. 

Not only does Um become negative for low tempera- 
tures, but also the net coupling coefficient 7 = 711 -t- 
722 ± 2712. Using Eqs. ([50| for detuned circular bands, 
we obtain for the s^ case: 



7+-=47rp(e^)T^ ^" /^ (51) 

Thus, 7+_ < for T* < 0.3^. Ahhough the sign of 

7-1 does not affect the criterion for phase coexistence, 

7^ < UmUs, it significantly changes the forms of the AFM 
and SC transition lines inside the coexistence region. In 
particular, a negative 7 implies that neither Tc nor M 
are suppressed in the AFM-SC coexistence regime - see, 
for instance, Eq. ([6]). 

A similar result for the AFM-SC coupling coefficient 
7 was obtained by Zhang et al. in the context of the 
cuprates^^. In a weak-coupling calculation at T = but 
finite disorder (otherwise the Matsubara sums would di- 
verge), they obtain a negative coupling coefficient be- 
tween a single-band d-wave SC order parameter and an 
itinerant AFM order parameter. Technically, the prob- 
lem of the competition between AFM and a single-band 
d-wave SC is equivalent to our two-band problem with 
tlie s^ SC state. Notice, however, that the coefficient 
7-1 only becomes negative at T* < T^^, i.e. the incom- 
mensurate AFM transition would happen before the cou- 
pling coefficient changes sign. Thus, 7-1 has no mean- 
ing in this regime and one would have to go back and 
calculate the coupling coefficient between an incommen- 
surate AFM order parameter and the SC order param- 
eter. However, the numerical calculations performed by 
Vorontsov et al. indicate that this coupling coefficient 
must be positiv o^^'^^ . Thus, in our approach, the most 
stable AFM state and superconductivity are always com- 
peting. 

Notice that this theoretical discussion about the in- 
commensurability of the AFM state does not take into 
account the coupling to the lattice degrees of freedom. 
As argued by many authors, an emergent nematic de- 
gree of freedom is present in the iron arsenides due to 
its frustrated magnetic structure^Sr— . The energy of the 
system is minimized by the onset of a nematic transition 
at Tncm ^ Tn. Due to the bilinear coupling between the 
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ncmatic order parameter and the shear distortion, the ne- 
matic transition is simuhaneous to a structural transition 
from the tetragonal to the orthorhombic phase. Key to 
this process is the commensurability of the magnetic fluc- 
tuations that give rise to this emergent nematic degree 
of freedoniS^. Thus, the inclusion of this extra degree of 
freedom could change the outcome of an incommensurate 
AFM state at low temperatures. 



D. Intraband pairing and Coulomb interaction 

In writing our weak-coupling expression for the SC in- 
teraction term, Eq. (|19p , we considered only an interband 
pairing interaction V = V12 = V21. If one includes ad- 
ditional intra-band pairing interactions Vn and V22, the 
only change in the free energy density, Eq. (|23p . is that 
the quadratic term —y (A*A2 -I- AjAi) is replaced by 
-J2ai3 i^~^)aB^a^li- This will of coursc change the 
gap equations, specially the value of Tc, and may also 
affect the ratio A1/A2. Yet, the inclusion of an intra- 
band pairing interaction will only change the quadratic 
Ginzburg-Landau coefficients a^^a, leaving the values of 
the quartic coefficients u^, u„i and 7 unchanged. Since 
our results regarding the coexistence or mutual exclusion 
between the SC and AFM states rely solely on the quartic 
coefficients, they will remain unchanged. 

Here, we assume that the interband pairing interac- 
tion V = V12 is originated from the coupling between 
electrons and collective modes of the system, such as 
phonons {V < 0) or paramagnons {V > 0), for exam- 
ple. With this in mind, we can investigate the effects of 
the electronic Coulomb repulsion by adding a renormal- 
ized Coulomb interaction U > 0. First, consider the case 
of an uniform Coulomb repulsion, with equal intraband 
and interband terms U . Formally, there is now a single 
interband interaction given by V + U, which is enhanced 
(reduced) in the case of s^ (s^^) pairing. Yet, due to 
the different origins of V and U, we here opt to write the 
total interband interaction in the form V + U. 

For the pure s"' state, it was previously shown that an 
uniform renormalized Coulomb interaction U is unable 
to completely destroy the SC state, i.e. Tc [U] never 
goes to zero, no matter the magnitude of the Coulomb 
interaction^. In order to demonstrate this, one writes 
the linearized gap equations in matrix form as Aq = 
Aq^A^ and analyzes the eigenvalues of: 
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Figure 11: Effective SC coupling constant Aeff = 
1/ In {W/aTa) as function of the ratio between the Coulomb 
repulsion and the pairing interaction U / \V\ for both a pure 
S+"'' state and a pure s^ state. Here, we considered the 
values |V| p\ = 0.1 and \V\ p2 = 0.3, but the conclusions are 
similar for arbitrary parameters. The inset is a zoom of the 
curve associated to the s^~ state. 



Tc, vanishes as C/ — ?> 1^^1/2, while in case of s^ -pairing 
the net pairing interaction stays finite even for infinite [/, 
where it holds A ([/ — > 00) = 2 \V\ P1P2/ (pi + 92)- These 
results are summarized in figure [TT] 

We next investigate the effects of an uniform Coulomb 
repulsion in the case where s^ -SC coexists with AFM. 
We now have: 



A = 



where 
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(53) 
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(54) 



with 



and the excitation energies are given by Eq. 
Ai = A2 — 0. The superconducting transition tempera- 
ture is again determined by the largest eigenvalue of A, 
which is given by: 



A 



-Upi -{V + U)p2 
-(!/ + [/) Pi -C/p2 



(52) 



The largest eigenvalue A determines the transition tem- 
perature through A^^ = In (W/aTc). For small U, it fol- 
lows that A = |y I ^p\P2 - \ (v^ ± 1/P2) ^ ■> where -1-1 
refers to s"*"^ pairing and —1 to s"' pairing, respectively. 
The suppression of the pairing interaction is significantly 
weaker for the s^ -state, in particular for similar den- 
sities of states p\ and p2- For s++-pairing A, and thus 



A = - (ri + r2 + 2ro) C/ - 2roy + (ri + ra + 2ro)^ C/ 



^2T^2 



-|-4rir2y^ + 4C/V (2rir2 -|- r^ri ^ ror2)] 



1/2 



(55) 



For C/ = 0, we find A = 2 {^/rir2 — rg) V. Let us 
assume that, in the absence of the Coulomb interaction, 
the system undergoes a SC transition at Tc- Imposing 
the vanishing of the largest eigenvalue, we obtain that 
A (Uc) = for Uc — —V/2, independent of the value of the 
magnetic order parameter M or of band structure details. 
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However, for the 5"' case, V > 0, implying Uc < 0. 
Therefore, the s^ SC state inside the antiferromagnetic 
phase is robust against an uniform Coulomb interaction, 
similarly to what happens for the nonmagnetic s"' state. 
We can also consider the case where the Coulomb in- 
teraction is not uniform, such that its value for intraband 
repulsion U is greater than its value for interband repul- 
sion U' < U. Now we need to determine the largest 
eigenvalue A of 



A 



{V + U') To - Ur2 ~{V + U') ri - Ur^ 
{V + U') r2 - Ura - {V + U') r^ - Un 



(56) 
A now vanishes for U — \V\ ±U' ^ where the plus (mi- 
nus) sign is to be used for the case of an s^ (■?'*"'') state. 
In this case, both states can be destroyed by a sufficiently 
large repulsion, yet the s++ state is destroyed easier than 
the s^ state. Note also that the condition U = \V\ztU' 
is the same for both situations of a pure SC state and a 
coexistent SC-AFM state. Thus, in general, magnetism 
does not seem to significantly influence the ability of the 
Coulomb repulsion to destroy the SC order. The renor- 
malization of the Coulomb interaction was also investi- 
gated using a renormalization group approach in Ref.^®. 
Even though the underlying reasoning is somewhat dif- 
ferent from our analysis, the conclusions of Ref.^^ are 
consistent with our result. The pair-breaking contribu- 
tion of the Coulomb interaction is less efficient for the 
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state, if compared to s^^-superconductivity. 



E. Reentrant Neel transition line and quantum 
fluctuations 



According to our phenomenological discussion in Sec- 
tion II, the strong suppression of the AFM order param- 
eter in the SC phase is also reflected in the reentrance 
of the AFM transition line, as shown by the calculated 
phase diagrams of figure [1] and confirmed by neutron 
diffraction measurements^*'. The same reentrant behav- 
ior is observed in some heavy fermions^^iM, where AFM- 
SC coexistence takes place as well. In the cuprates, the- 
oretical models also proposed that a similar reentrance is 
present in the phase diagram—.. 

Let us investigate in more detail the form of the reen- 
trance line in the iron arsenides. For simplicity, we follow 
Ref.''^ and consider a small perturbation of the particle- 
hole symmetric band structure, ^2,k+Q = ~^i,k ~ '^S^, 
with S^ — Sq + 82 cos 2iy9, as we explained in Section IV- 
B. Assuming that Sq and 82 satisfy the conditions for 
coexistence between AFM and s^ -SC, we expand the 
free energy (P5)) only in powers of M, keeping the SC 
gap A = Ai = — A2 fixed. This last assumption is justi- 
fied at low temperatures, where the SC order parameter 
saturates (see Fig. |9]). For T ^ 5 <C A, we obtain for 
the quadratic magnetic coefficient: 



a™(A,r) 



A -A, 



^^1^ -* <"' 



with (5-^ = ((5^) — &\ + <^2/2- Here, Ac corresponds to 
the T = value of the SC gap where the quantum phase 
transition from the superconducting normal state to the 
superconducting state with antiferromagnetic long range 
order takes place. Eq. ([57)1 implies that, at T = 0, there 
is AFM order for A < Ag. The negative sign in front of 
the temperature dependent term also implies that Tf^ is 
finite for A > Ac. Therefore, it correctly captures the 
reentrance of the AFM line. 

The presence of the exponential term e"T is a conse- 
quence of the fact that, inside the SC state, quasiparticle 
excitations are fully gapped. Due to this exponential de- 
pendency, the reentrant Tjq line approaches the quantum 
critical point with an exponentially steep slope, i.e. as 
an almost vertical line. This is in agreement with our 
calculated phase diagrams from figure [TJ as well as with 
the phase diagrams obtained by Vorontsov et al^. 

These results were derived using a mean-field ap- 
proach. Close to T = 0, the presence of quantum fluc- 
tuations change this scenario. To illustrate their effects, 
we consider collective magnetic fluctuations in the vicin- 
ity of this quantum critical point. In order to properly 
describe long range mangnetic order we have to include 
interlayer coupling and consider an effective three dimen- 
sional quantum rotor model. Due to the fact that quasi- 
particle excitations are gapped, we consider a rotor model 
with dynamic critical exponent^ z = 1. Within a self 
consistent large- A^ theory, where N refers to the number 
of components of the rotor, we obtain a renormalization 
of the coefficient in am — > Om due to critical fluctiations: 



-^m ^m 
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dh 



{2nf fir; 



(58) 



An expansion at low temperatures yields cim — o,m = 
CT^, with C > 0, which dominates over the exponential 
term e~T that follows from the mean field theory. The 
presence of power-law corrections in am is more general 
than our self-consistent large- A^ theory, and is expected 
to occur due to the presence of massless critical fluctua- 
tions. Due to the fact that the coefficient C is positive, 
such quantum fluctuations suppress the magnetic reen- 
trant behavior at very low temperatures. However, as 
usual in systems in the weak coupling regime, the crit- 
ical region where quantum fluctuations are relevant is 
expected to be very small, and probably hard to be de- 
tected experimentally. 
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V. LOCALIZED VERSUS ITINERANT 
MAGNETISM 

A key conclusion of our calculation is that homoge- 
neous coexistence of superconductivity and magnetism 
is only allowed in case of unconventional s^ -pairing 
state, while both ordered states exclude each other in 
case of conventional s++-pairing. This conclusion seems 
to be at odds with the well known fact that antiferro- 
magnetism and superconductivity do coexist homoge- 
neously in a number of materials where the evidence 
for conventional electron-phonon pairing is very strong, 
such as the borocarbides^^ RNi2B2C and the ternary 
superconductors^^ RMogSg and RRh4B4, with R denot- 
ing a rare earth. The crucial difference between these rare 
earth based systems and the iron pnictide superconduc- 
tors is that the magnetism in the former is due to local- 
ized rare earth spins while in the latter the same electrons 
that superconduct are responsible for the entire ordered 
moment. Thus, for our argumentation in the pnictides to 
hold, it is essential that the same electrons that form the 
Cooper pair condensate are responsible for the ordered 
moment. This is evident from our Hamiltonian, Ea. (|lip . 
where the order parameters A and M are expectation 
values of electronic states of the same bands. This is the 
reason for the highly symmetric form of the free energy 
of Eas. pil32p and why the Ginzburg-Landau coefficients 
for the quartic magnetic, superconducting and coupling 
terms are closely related to each other. 

In order to demonstrate explicitly that the phase dia- 
gram of competing magnetism and superconductivity is 
very different in case of localized spins, here we analyze 
this problem in some detail. We recall that the total 
Hamiltonian is given by 'Hq+'Hafm+'Hsc- We keep the 
same terms for the kinetic and superconducting parts, 
given by Eqs. P^ and (|19p . respectively. The pairing 
interaction V in "Hsc niight, for example, be due to the 
electron-phonon interactions. At this point we are not 
concerned whether systems like the RNi2B2C are indeed 
characterized by a corresponding two band model. In- 
stead, we are primarily interested in comparing localized 
and itinerant magnetism for a system with otherwise un- 
changed electronic structure. It will become evident be- 
low that our analysis is in fact more general. The crucial 
new term in the Hamiltonian is Hafm which is replaced 
by: 



than the corresponding Kondo temperature. To proceed, 
we perform a mean field analysis of this model. We in- 
troduce the expectation values: 
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(60) 



with magnetic ordering vector Q. For definiteness, we 
consider Jk > 0, implying that (S'f ) and (sf) have op- 
posite sign. Since we ignore the Kondo effect, our final 
results are independent on the sign of Jk- 

In analogy to the theory of itinerant magnetism we 
perform a mean field calculation, giving rise to the total 
free energy density F = Fs+Fo,c\+Fsc with contributions 
from localized spins, Fs, from the SC condensate Fsc, 
and from the electronic part, i^o.ci, which also includes 
the order parameters coupling. The last two terms are 
completely analogous to the case of an itinerant AFM 
state competing with SC, Eq. (P5)) . if we identify the 
magnetic order parameter as: 
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JKmioc 



(61) 



Recall that M in our notation is the antiferromagnetic 
potential that causes a gap for Bragg reflected points of 
the Fermi surface. In case of itinerant magnetism this 
gap is due to the electron-electron interaction / and the 
moment of the itinerant electrons. Now the microscopic 
origin of M is very different. Yet, the expression for the 
energy of the conduction electrons is still given by: 
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(62) 



with the same excitation energies Ea,k = 
Ea,kiAi,A2,M) from Eq. ^. The contribution 
to the energy due to the pairing interaction is unchanged 
as well and given by 



Fsc = -J2^-;aia^ 



(63) 
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(59) 



Finally, the free energy density due to the localized 
spins is: 



Here Si refers to a localized spin-S* operator and Jk is 
the exchange coupling between localized spins and con- 
duction electrons. We are interested in the regime where 
Jk leads to magnetic long range order via the RKKY 
mechanism with Jrkky (r) — JkXs (r), where Xs (r) de- 
notes the electronic spin susceptibility. In the regime of 
antiferromagnetism with large ordered local moments, it 
is possible to neglect the Kondo effect as Jrkky is larger 



s 
Fs = -T\n Y^ e™^'' + Jifmiocmei, (64) 

where h = Jkitt-cI is the Weiss field of a single spin-S". 
Since both magnetizations mioc and TOci order simulta- 
neously, we can eliminate rrioi and express the Landau 
expansion in terms of M. To this end we use ttiioc = 
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- dFs,/dh\^^j^^^^ and solve for TOoi = md (mjoc). Using 
M of EqinU instead of mioc , we find to leading order 
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After inserting this result for rrici into the free en- 
ergy, Eq |64[ we can expand it and thus determine the 
Ginzburg-Landau expansion simultaneously for the SC 
and AFM order parameters, Ai, A2 and M. We obtain 
the exact same expressions for the coefficients related to 
the superconducting order parameter, a^^ap and Ms.a, as 
well as for the coupling 7q,^ between the AFM and SC or- 
der parameters. Despite the same formal expression, the 
physical interpretation of the order parameter coupling 
terms 7^^ is somewhat different now. It reflects changes 
in the conduction-electron mediated RKKY interaction 
due to the onset of superconductivity. 

The only difference in the values of the Ginzburg- 
Landau parameters due to the presence of localized spins 
is for the coefficients of the magnetic order parameter. 
We find: 
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where u^ is the quartic coupling of itinerant spins given 
by Ea ((35|l . The additional term in Eq. ([68| is solely de- 
termined by the Neel temperature. Tat, the size of the 
spin, iS*, and the coupling Jk- TV is determined via 
flm {Tm) = and given, as expected, by the RKKY cou- 
pling 



T^w = - asficXph (Q) 



(69) 



with bare spin susceptibility of the conduction electrons 
at the ordering vector Q. Since Eq. (gS]) gives u^ w 
p/ (2T^) , it follows that the relative change in the mag- 
netic Ginzburg-Landau coefficient is: 



^(JKp)^ln3(a,V^,V^). (70) 

as 



The additional logarithmic term In^ (a^^ J^^p~^) oc- 
curs only near particle-hole symmetry and is replaced 
by a constant of order unity away from particle hole- 
symmetry. The prefactor Ps las grows as S'^ for large S. 
Thus, it is easily possible that the quartic coefficient of 



the magnetic order parameter is significantly enhanced 
in case of localized spins. For example, the relative cor- 
rections are around 200% for S = 7/2 and JkP — 0.025. 
Since the order parameter coupling and the quartic coeffi- 
cients of the superconducting term are unchanged, it fol- 
lows that the condition 7 < ^/UsUm for coexisting order 
can now be fulfilled easier than in the case of purely itin- 
erant systems. This offers a natural explanation for the 
observation of homogeneous coexistence of both phases 
in systems such as the RNi2B2C and addresses the fact 
that coexistence observed in systems with localized spins 
is not in contradiction to our conclusions. For complete- 
ness, we also analyzed a model with additional magnetic 
interactions between localized spins that are not captured 
by the RKKY mechanism, adding to "Hafm of Ea. ([CT|) 
the term ^ X^i Jij^i ' Sj. This new term will change the 
value of Tn , but not affect the expression Eq. (|68|) for Um ■ 

Additional consequences for localized spins are that 
the coefficient am,o = {asJK) of the temperature de- 
pendent quadratic coefficient «„ = am^o {T — Tn) is ex- 
pected to be larger compared to the corresponding coef- 
ficient Os^o — p/Tc of the superconducting order param- 
eter if we consider the multicritical point Tc — T/v near 
particle-hole symmetry, since amfi/as,o — Xph (Q) I P- 
With Eqini follows then that it becomes harder to achieve 
a suppression of the magnetization with d'M.'^ /dT > 
below the superconducting transition. The observed 
suppressioE^i^ii^i^ of M in the coexistence region of 
Ba (Fei_xCox)2 As2 and Ba (Fei_xRhx)2 As2 is therefore 
an indication that the same electrons are responsible for 
both states and that magnetism is itinerant in these ma- 
terials. On the other hand the condition Qs^ol < o-mfiUm 
for suppression of T^ in the magnetically ordered state 
can easily be fulfilled. Thus, while SC in systems with 
localized spins is affected by magnetic long range order, 
the opposite does not seem to hold and AFM is rather 
indifferent to SC. 

Finally we comment on the relevance of this calcula- 
tion for heavy fermion superconductors, such as CeRhIn5 
and UPt3. In CeRhlns, the coexistence between mag- 
netism and superconductivity has been investigated in 
great detaiiSS,, while in UPta there is clear evidence^i for 
suppresion of magnetism below T^. The heavy fermion 
system are believed to be properly described by the 
Kondo lattice Hamiltonian^ with coupling between lo- 
calized and conduction electrons as in Hafm of Eq. (|59p . 
However, our analysis of this model, where we completely 
ignored the Kondo effect and the emergence of a heavy 
electron state is inadequate for such systems. In fact 
one expects that a system in the heavy electron state 
is better described by the theory employed here for the 
FeAs systems, yet the interactions / and V as well as the 
quasi-particles masses are heavily renormalized due to 
Kondo lattice screening. Thus, while a detailed theory 
for the competition of magnetism and superconductiv- 
ity in heavy electron states is complex, we do expect a 
similar competition in the FeAs systems and these heavy 
fermion compounds. The very similar behavior of the SC 
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and AFM transition lines in the phase diagrams^SiS^iS^ of 
Ba (Fei_xC0x)2 As2 and CeRhlns, as well as the suppres- 
sion of the magnetizatioE^ of UPta below Tc, certainly 
support this view. 



VI. SUPPRESSION OF SC IN THE 
OVERDOPED REGION 

So far we have analyzed the competition between AFM 
and SC with the consequent suppression of the super- 
conducting state in the underdoped region of the phase 
diagram of the iron arsenides. In the overdoped region 
there is no magnetically ordered state, yet SC is also sup- 
pressed and eventually disappears. The Fermi surface in 
this part of the phase diagram also changes significantly: 
for electron (hole) doped samples, it is characterized by 
increasingly large electron (hole) pockets and decreas- 
ingly small hole (electron) pockets^'^, which eventually 
disappear at a certain doping level. In this section, we 
investigate how the the disappearance of these pockets 
from the Fermi surface affects the transition temperature 
of the pure SC state. 

We consider, once more, one hole pocket centered at 
the Brillouin zone and one electron pocket displaced by 
Q from the zone center. For definiteness, we use the 
band dispersions (J13p with Eq = ei^o = £2,0 ^-nd vary the 
chemical potential /i. Here, we consider the effects of 
electron doping only, such that fj, > 0; the case of hole 
doping (/z < 0) is analogous and the same conclusions 
hold. For simplicity, we first consider V to be constant 
as /i changes and neglect intraband pairing interactions. 
The linearized gap equations are then given by: 



Ai = - 



VA2 



rW 



A, 



J-w 
2 J-w 






tanh 



V^Tej 



tanh 



V^Tej 



e 



(71) 



where, once again, W denotes an upper energy cutoff as- 
sociated to the pairing interaction. In a two-dimensional 
system, the density of states pi (^) is constant if the en- 
ergy f falls inside the bands. Therefore, diagonalization 
of the linearized gap equations give the following implicit 
expression for Tc'. 



A§ 



dx 
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2Tc 
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, tanh a: tanhu ,„„, 

dy ^ (72) 

fo+ii") -ixy 



where we introduced Ag = Vy/pip2- For simplicity, let 
us first consider the special case W = So', the main con- 
clusions hold for an arbitrary W. 

Although a complete analytical solution for Tc (/i) is 
not available from Eq. ([7^ . we can obtain some impor- 
tant limits. For small /x ^ Eq, we obtain that T^ decreases 
linearly with respect to Tc = Tc (fi = 0): 



Tc = Tjo) f 1 - -^ 

4£o 



(73) 



Equivalently, we can show that the effective coupling 
( T" ^%~ ) decreases linearly with p. 
At the special point p — Sq, 



constant Ap 



In 



from its /i = value A. 
i.e. when the hole pocket shrinks to a single point in the 
Fermi surface, the effective coupling constant is reduced 
to about 70% of its initial value, Acff = A/-\/2, implying: 



Tc = r(o)e 



V^-1 



(74) 



Even though the vanishing of the hole pocket does not 
cause Tc to vanish, it does signal the onset of a regime 
where the coupling constant decays very strongly with 
respect to the chemical potential p, in contrast to the 
case of small chemical potential, where the decay Acff was 
linear in p. To illustrate this point, consider the regime 
of eo < A* < 2eo such that p — Eq ^ Tc. Notice that this 
condition is not too restrictive, since Tc < Tc <$^ Ea, by 
construction. Then, it follows that: 



Aeff = -:- In 



£0 



p- Eo 



(75) 



implying: 



Tc = ri°) 



exp 



A2ln 



M-eo 



(76) 



Clearly, the effective coupling constant only vanishes 
a.t p = W + Eq = 2eq, i.e. where the range of energies for 
which the net attractive interaction is positive does not 
cross the hole pocket. However, it is already significantly 
reduced for values of the chemical potential much smaller 
than that one. In figure[T2l we present the behavior of Tc 
for two values of the effective zero-doping coupling con- 
stant, Ao = 0.2 and Aq = 0.3. Notice that, in both cases, 
Tc decreases moderately for p < Eq (i.e. when the hole 
pocket is small but still present) and then is strongly re- 
duced for p > Eq (i.e. when the hole pocket disappears 
from the Fermi surface) . The same conclusions also hold 
for the case of an arbitrary cutoff W > e^: the special 
point of the phase diagram where the hole pocket is re- 
duced to a Fermi point marks the onset of a dramatic 
reduction of Tc, no matter the initial value of the cou- 
pling constant Aq. 

As we mentioned, Tc actually only vanishes at /i = 
W -\- E{). However, in this region where it is strongly 
suppressed, even an uniform Coulomb interaction is able 
to completely destroy superconductivity. This is to be 
contrasted to the optimally doped region, where the s^ 
state is robust against an uniform Coulomb repulsion, as 
we discussed in the previous section. 



19 



T /T(0) 
C ' C 




X /j(0) 

n ' r 



0.2 0.4 0.6 0.8 1.0 |i/8„ 



c ' c 




0.2 0.4 0.6 0.8 1.0 |^/8o 



Figure 12: Superconducting transition temperature Tc (in 
units of its value at zero chemical potential Tc ) as function 
of the chemical potential /i (in units of £o) for both an initial 
effective SC coupling Ao = 0.2 (a) and Ao = 0.3 (b). The 
insets show the changes in the Fermi surface with fi, with the 
red circle (blue ellipse) denoting the hole (electron) pocket. 
In the upper right corner of each panel we present a zoom of 
the region around fi — eo where the electron pocket is reduced 
to a Fermi point. 



To illustrate how the s^ SC state is killed in the re- 
gion of strong suppression, /i — eo ^ Tc, consider W — Eq 
and an uniform repulsion U > 0. For simplicity, we fo- 
cus only on the limit oi U ^ V. To leading order, the 
equation determining the effective SC coupling constant 
is given by: 
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Defining Xi = Vpi, we obtain: 
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Note that Acff vanishes at /i* = eg (l + e ^^^^), but is 
positive for fj, < fi*. Tlius, a sufficiently large Coulomb 
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Figure 13: Comparison between the behavior of Tc (in units 
of its value at zero chemical potential Tc ) as function of the 
chemical potential fi (in units of eo) for a two-dimensional 
and a three-dimensional system (green and red curves, re- 
spectively) . The coupling contants were chosen to give Tc ~ 
0.03 W in both cases, with W — eo- 



contrary to what happened when the hole pocket was 
part of the Fermi surface, where AoS — > const, even as 

U ^)- 00. 

The previous analysis holds for a 2D system. In 
the case of a three-dimensional system, the density of 
states is not constant anymore, but given by pi (^) = 
/I — ^ for the hole band and by p2 {£,) — 
^ for the electron band. Here, we neglected 



■/i- 



interaction U is now able to destroy the 



SC state. 



ci\/io~ 
C2VS0' 

the anisotropy of the electron band and denoted unim- 
portant constants by Ci . The main difference from the 
two-dimensional case is that the densities of states go to 
zero at the bands edges, leading to a significant stronger 
reduction of Tc as the chemical potential increases. This 
is illustrated in figure I13[ where we compare the solu- 
tions of the linearized gap equations ([7T|) for the cases of 
a 2D and a 3-D system, with coupling constants chosen 
to yield Tc (^ = 0) « 0.03 W. Thefore, close to the point 
of the phase diagram where the hole pocket disappears, 
the space dimensionality matters much more than in the 
usual Cooper problem. 

Our results suggest that the main factor responsi- 
ble for the complete suppression of the SC state in 
the overdoped region of the iron arsenides phase dia- 
gram is the evolution of the Fermi surface with doping. 
This is corroborated by ARPES measurements-^"'"'^ on 
Ba (Fei_xCox)2 As2, winch show no superconductivity in 
tlie overdoped region after the Fermi surface loses one of 
its pockets around x ~ 15%. The fact that superconduc- 
tivity requires the presence of the hole pocket and the 
observation that the superconducting gap on the hole 
and electron pockets are very similar, strongly support 
the view that the pairing interaction in the pnictides is 
due to interband coupling. 
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VII. CONCLUDING REMARKS 

We showed that the abihty of superconductivity and 
antiferromagnetism to order simultaneously depends sen- 
sitively on the nature of the Cooper pair wave function. 
In a two-band system with particle-hole symmetry, the 
itinerant AFM state and the unconventional s^ SC 
state are exactly in a borderline regime between phase 
coexistence and phase separation. In contrast, the con- 
ventional 5+"^ state is deep in the regime of mutual phase 
exclusion. We further demonstrated that this result holds 
regardless of additional details of the band structure or 
the system's dimensionality. It does not change either 
when one considers the presence of intraband pairing 
interactions or the effects of Coulomb repulsion. The 
robustness of this result, valid around the multicritical 
point T/v — Tc, is related to the quasiparticle excitation 
spectrum of the system, which depends on the peculiar 
combination M^ + |A| = N^ only. This is the root of 
the SO (5) symmetry of the free energy expansion, which 
has been shown to hold not only in the mean-field level, 
but also in the strong coupling limit, as a subgroup of an 
emergent SO (6) symmetry present in the Hamiltoniait^. 
Furthermore, the inclusion of fluctuations in our free en- 
ergy expansion are known to not change the condition 
for having the borderline regimej^i*^ g — 0. All these 
facts suggest that this simple result is a much more gen- 
eral property of this type of system. We also demon- 
strated that for our results to hold, it is crucial that the 
same electrons that form Cooper pairs are responsible for 
the formation of the ordered moment. For instance, we 
showed that an AFM state generated by localized mo- 
ments has a tendency of being indifferent to SC, falling 
much easier in the regime of phase coexistence. 

When applying this result to the pnictides, one has 
to critically evaluate whether the assumptions that were 
made are too restrictive. We assumed a certain vicinity 
to particle-hole symmetry. It is evident that the real ma- 
terials do not have perfectly nested bands^'^. Yet, there 
are clearly two sets of pockets whose centers are separated 
by the ordering vector Q. As we showed, small pertur- 
bations in the ellipticity of one pocket or in the chemical 
potential bring the s"' state to the regime of mutual 
exclusion. However, simultaneous perturbations in both 
quantities can bring the system to either regime, as our 
figure [5] illustrates. Most importantly, small deviations 
from perfect nesting are not sufficient to lead to simulta- 
neous order of magnetism and s''"+-pairing close to the 
multicritical point. The fact that particular details of the 
band structure are able to bring the s"* state either to 
the phase coexistence or to the phase separation regime 
is, in our view, what makes some of the iron arsenides dis- 
play second-order AFM-SC transition and others, first- 
order AFM-SC transition. For Ba (Fei_xCox)2 As2, we 
were able to independently "fit" the effective band struc- 
ture parameters to give the correct doping and Tj\j de- 
pendence of the zero-temperature magnetization in the 
absence of superconductivity^". Having this set of pa- 



rameters, which did not depend on any SC property, 
we were able to obtain the complete mean-field x — T 
phase diagram and to readily calculate the parameter 
g^ verifying that, as evidenced by many experimental 
probeS )^'^'^'^" — Ba (Fei_xCox)2 As2 has homogeneous co- 
existence between SC and AFM states. 

For other compounds that are believed to dis- 
play phase separation, like LaFeAs (Oi_xFx) and 
(Bai_xKx) Fe2As2, we do not have the same systematic 
diffraction measurements that allow a reliable extraction 
of effective two-band parameters. Even though tight- 
binding fits to DFT-calculated band structures are avail- 
able, they usually refer to the parent compounds. The 
problem to extrapolate them to finite doping is that, in 
these materials, doping is not on the Fe site nor on the 
FeAs plane, which makes it much more difficult to make a 
direct association between electronic occupation number 
and doping level, as we did for Ba (Fei_xCox)2 As2. Fur- 
thermore, in some of the pnictides, the issue of whether 
they have homogeneous or heterogeneous AFM-SC co- 
existence is still not completely settled. For exam- 
ple, in (Bai_xKx) Fe2As2, while many experiments find 
phase separation, Mossbauer spectroscopy identifies mi- 
croscopic coexistence^*'. In SmFeAs (Oi_xFx), muon spin 
rotation® finds homogeneous coexistence while a work 
combining powder x-ray diffraction, Mossbauer spec- 
troscopy and nuclear resonant forward scattering finds 
phase separation^!. 

In what concerns the pairing state, although many the- 
ories for pairing due to electronic interactions predict a 
nodeless s^ state, other models also suggest that acci- 
dental (i.e. not related to symmetry) node o^^i^^i^® or even 
d-wave nodal states^ could be present. Furthermore, 
penetration depth experiments are in principle consistent 
with the existence of accidental nodes^- . The inclusion 
of gap nodes in our model, as we discussed, tend to move 
the SC state from the borderline regime to the mutual 
exclusion regime, although not as deep as in the s^^ 
case. Clearly, Aig sign-changing gap functions with spe- 
cific configurations of nodes could be as effective as small 
perturbations of the particle-hole symmetric band struc- 
ture in making the AFM-SC state go from one regime to 
the other. Thus, while we cannot discard the presence 
of unconventional nodal states coexisting with AFM, we 
can certainly discard the conventional s++ state: it is 
simply incompatible with magnetism. 

Regarding the nature of the magnetic state, there is 
still some debate whether the AFM phase is due to con- 
duction electrons or localized spins. Our analysis shed 
light on this subject: as our calculations demonstrated, 
a magnetically ordered phase with localized spins can 
coexist much easier with superconductivity. This holds 
even for conventional BCS states, like in the case of the 
ternary^^ and quaternary^ rare earth compounds. Fur- 
thermore, the magnetization is not so affected by the SC 
condensate as in the case of purely itinerant magnetism. 
Therefore, the experimental observation of reentrance of 
the paramagnetic phase inside the SC dome^S seems to 
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rule out the AFM state formed only by localized mo- 
ments. 

As we showed, in the underdoped side of the FeAs 
phase diagram, both T/y and Tc are suppressed due to 
the competition between the AFM and SC phases. The 
Neel transition line is even bent back, approaching the 
X axis vertically, similar to the case of a first-order tran- 
sition line. The suppression of Tc is milder; this differ- 
ence is probably due to the fact that the AFM gap opens 
only around the corresponding Bragg scattered points,— 
whereas the SC gap opens isotropically around the Fermi 
surface. In the overdoped side, however, SC is the only 
thermodynamic ordered phase of the system and yet it 
is suppressed similarly to the underdoped side. Our cal- 
culations using the coupled SC gap equations shows that 
this suppression is related to the disappearance of one 
of the pockets from the Fermi surface. As suggested by 
our figures [1^ and [T31 and by ARPES measurements^"— 
in Ba (Fei_xCox)2 As2, the main cause of suppression of 
SC in the overdoped side seems to be the doping-induced 
changes in the Fermi surface rather than possible changes 
in the magnitude of the pairing interaction. Interestingly, 
the Fermi surface is also indirectly responsible for the 
suppression of Tc in the underdoped side, since it is the 
driving force of the AFM instability'^^. 

The calculations we presented here did not take into ac- 
count another important degree of freedom present in the 
iron arsenides, the orthorhombic distortion. However, in 
separate works, we showed that the competition between 
SC and AFM, combined with the coupling between ne- 
matic degrees of freedom and structural distortion, are 
able to consistently explain the observed back-bending 
of the structural transition line^'^ inside the SC dome as 
well as the increase of the shear modulus®^ below Tc- 

Finally, we comment that in many models of uncon- 
ventional pairing in the iron arsenides, the bosons re- 
sponsible for the formation of Cooper pairs are spin 
fluctuations^"—!^"—. In this case, the pairing interac- 
tion V itself would also be sensitive to the presence of 
magnetic long range order. Clearly, this would change 
some details of the coexistence state, particularly the 
form of the transition lines. Yet, our main conclusions 
still hold in this case, since the decision about the coexis- 
tence between AFM and SC is made when both transition 
temperatures are close, implying that the order parame- 
ters are only infinitesimal. 



Appendix A: Diagrammatic interpretation of the 
Ginzburg-Landau coefficients 

Here, we rederive the Ginzburg-Landau expansion in 
Eq. (|23p by explicitly integrating out the fermionic de- 
grees of freedom of a system with competing SC and 
AFM. This method is useful since it provides a diagram- 
matic interpretation for the coefficients, as presented in 
figure 131 We first define the Nambu operator: 
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Ckt C_ki «k+Qt "-k-Qi 
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and generalize the uniform order parameters Aq, and M 
to inhomogeneous functions of space and time A^j (k.w,,) 
and M(j.^^), where a;„ — (2n -|- l)7rT is a fermionic Mat- 
subara frequency. Thus, denoting k — (k, Wn), we obtain 
the Green's function: 



^kX' = (*^" - ?k j 5k.k' - Uk-k' (A2) 

where the hat denotes a matrix in Nambu space and: 



Uk-k' = ^k-k' + Mk-k' 



(A3) 



With the help of the Pauli matrices r, we can write 
the 4x4 Nambu matrices as 
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as well as 
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Therefore, after introducing the condensation energy 
of the magnetic and superconducting phases, we obtain 
the action: 
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where J^, = TJ2J2 ^i^d /^ = h ! d'^^ !o dr, with v 

n k 

denoting the volume of the system. The fermions can now 
be integrated out, yielding an effective action in terms of 
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the collective AFM and SC fields S'off. The partition 
function of the free fermions is given by: 



Zo = det (~Qq ^ 



(A8) 



where we defined the non-interacting Nambu Green's 
function Gq^-i^i = (iw„ — ^k ) <^fe.fe'- The effective action 



reads: 
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Here, the trace refers to the sum over momentum, fre- 
quency and Nambu indices. Notice that, for uniform 
AFM and SC gaps, the total free energy density (free- 
fermions contribution included) is given by: 
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Since det ( G^ M — J2 i^n + ^a it.) ^ ^^ '^^'^ evaluate 



the Matsubara sum to obtain: 
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(All) 
with excitation energies Ea.k given by the positive roots 
of (221) ■ Minimization of the free energy (JA11[) with re- 
spect to Aq, and M leads then to the gap equations (I^Ul) . 
Going back to the effective action (IA9[) . we can now 
perform a Ginzburg-Landau expansion in the AFM and 
SC order parameters. Expansion of the logarithm yields, 
for the free energy relative to the paramagnetic, normal 
phase: 
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Considering static and homogeneous AFM and SC 
gaps and performing the traces in the Nambu space, we 
obtain the free energy density expansion of Eq. (P5)) with 
Ginzburg-Landau coefficients: 
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In the previous expressions, we introduced the non- 
interacting single-particle Green's functions Gi.fe — 

{iujn -^i,k)"^ and G2,k = (iw„ - ^2,k+Q)~^ Thus, the 
coefficients depend only on the band structure and on the 
magnitude of the electronic interactions. The Feynman 
diagrams associated to the quartic coefficients of (IA13P 
are presented in figure 21 Evaluation of the Matsubara 
sums then leads to ([55l) . ([571) and (|551) for the Ginzburg- 
Landau coefficients. 



^ Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 
Am. Chem. Soc. 130, 3296 (2008). 

2 X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, and D. 
F. Fang, Nature 453, 761 (2008). 

3 G. F. Chen, Z. Li, D. Wu, G. Li, W. Z. Hu, J. Dong, P. 
Zheng, J. L. Luo, and N. L. Wang, Phys. Rev. Lett. 100, 
247002 (2008). 

'^ M. Rotter, M. Tegel, and D. Johrendt, Phys. Rev. Lett. 
101, 107006 (2008). 

^ H. Luetkens, H.-H. Klauss, M. Kraken, F. J. Litterst, T. 
Dellmann, R. Klingeler, C. Hess, R. Khasanov, A. Am- 
ato, C. Baines, M. Kosmala, O. J. Schumann, M. Braden, 
J. Hamann-Borrero, N. Leps, A. Kondrat, G. Behr, J. 
Werner, and B. Bi;,oechner, Nature Mater. 8, 305 (2009). 

^ A. J. Drew, Ch. Niedermayer, P. J. Baker, F. L. Pratt, S. 
J. Blundell, T. Lancaster, R. H. Liu, G. Wu, X. H. Chen, 
I. Watanabe, V. K. Malik, A. Dubroka, M. Ri^oessle, K. 



W. Kim, C. Baines, and C. Bernhard, Nature Mater. 8, 
310 (2009). 

N. Ni, M. E. Tillman, J.-Q. Yan, A. Kracher, S. T. Han- 
nahs, S. L. Bud'ko, and P. C. Canfield, Phys. Rev. B 78, 
214515 (2008). 

J.-H. Chu, J. G. Analytis, C. Kucharczyk, and I. R. Fisher, 
Phys. Rev. B 79, 014506 (2009). 
C. R. Rotundu, D. T. Keane, B. Freelon, S. D. Wilson, 

A. Kim, P. N. Valdivia, E. Bourret-Courchesne, and R. J. 
Bhgeneau. .arXiv: 0907.1308 (2009). 

T. Goko, A. A. Aczel, E. Baggio-Saitovitch, S. L. Bud'ko, 
P. C. Canfield, J. P. Carlo, G. F. Chen, P. Dai, A. C. 
Hamann, W. Z. Hu, H. Kageyama, G. M. Luke, J. L. Luo, 

B. Nachumi, N. Ni, D. Reznik, D. R. Sanchez-Candela, A. 
T. Savici, K. J. Sikes, N. L. Wang, C. R. Wiebe, T. J. 
Williams, T. Yamamoto, W. Yu, and Y. J. Uemura, Phys. 
Rev. B 80, 024508 (2009). 



23 



^^ H. Fukazawa, T. Yamazaki, K. Kondo, Y. Kohori, N. 

Takeshita, P. M. Shirage, K. Kihou, K. Miyazawa, H. Kito, 

H. Eisaki, and A. lyo, J. Phys. Soc. Jpn. 78, 033704 (2009). 
^2 J. T. Park, D. S. Inosov, C. Niedermayer, G. L. Sun, D. 

Haug, N. B. Christensen, R. Dinnebier, A. V. Boris, A. 

J. Drew, L. Schulz, T. Shapoval, U. Wolff, V. Neu, X. P. 

Yang, C. T. Lin, B. Keimer, and V. Hinkov, Phys. Rev. 

Lett. 102, 117006 (2009). 
^^ Y. Laplace, J. Bobroff, F. Rullier-Albenque, D. Colson, 

and A. Forget, Phys. Rev. B 80, 140501(R) (2009). 
^* M. -H. Julien, H. Mayaffre, M. Horvatic, C. Berthier, X. 

D. Zhang, W. Wu, G. F. Chen, N. L. Wang, and J. L. Luo, 

Europhys. Lett. 87 37001 (2009). 
1^ C. Bernhard, A. J. Drew, L. Schulz, V. K. Malik, M. 

Rossle, C. Niedermayer, T. Wolf, G. D. Varma, G. Mu, 

H. H. Wen, H. Liu, G. Wu, and X. H. Chen, New J. Phys. 

11, 055050 (2009). 
" F. Ning, K. Ahilan, T. Imai, A. S. Sefat, R. Jin, M. A. 

McGuire, B. C. Sales, and D. Mandrus, J. Phys. Soc. Jpn. 

78, 013711 (2009). 
^^ C. Lester, J.-H. Chu, J. G. Analytis, S. C. Capelh, A. S. 

Erickson, C. L. Condron, M. F. Toney, L R. Fisher, and S. 

M. Hayden, Phys. Rev. B 79, 144523 (2009). 
^^ D. K. Pratt, W. Tian, A. Kreyssig, J. L. Zarestky, S. 

Nandi, N. Ni, S. L. Bud'Ko, P. C. Canfield, A. I. Gold- 
man, and R. J. McQueeney, Phys. Rev. Lett. 103, 087001 

(2009). 
^^ A. D. Christianson, M. D. Lumsden, S. E. Nagler, G. 

J. MacDougall, M. A. McGuire, A. S. Sefat, R. Jin, B. 

C. Sales, and D. Mandrus, Phys. Rev. Lett. 103, 087002 

(2009). 
^° R. M. Fernandes, D. K. Pratt, W. Tian, J. L. Zarestky, 

A. Kreyssig, S. Nandi, M. G. Kim, A. Thaler, N. Ni, P. 

C. Canfield, R. J. McQueeney, J. Schmalian, and A. L 

Goldman, Phys. Rev. B 81, 140501(R) (2010). 
2^ W. Z. Hu, J. Dong, G. Li, Z. Li, P. Zheng, G. F. Chen, 

J. L. Luo, and N. L. Wang, Phys. Rev. Lett. 101, 257005 

(2008). 

^2 R. M. Fernandes and J. Schmahan, larXiv:1005.2174l 

(2010). 
^^ D. J. Singh and M.-H. Du, Phys. Rev. Lett. 100, 237003 

(2008). 
^^ V. Cvetkovic and Z. Tesanovic, EPL 85, 37002 (2008). 
^^ M. M. Korshunov and I. Eremin, Phys. Rev. B 78, 

140509(R) (2008). 
^^ P. M. R. Brydon and C. Timm, Phys. Rev. B 80, 

174401(2009). 
^^ I. Eremin and A. V. Chubukov, Phys. Rev. B 81, 024511 

(2010). 
^^ W. Baltensperger and S. Strassler, Phys. Kondens. Mater. 

1, 20 (1963). 
^^ L. N. Bulaevskii, A. I. Rusinov and M. Kulic, J. Low Temp. 

Phys. 39, 255 (1980) 
^° M. J. Nass, K. Levin, and G. S. Crest, Phys. Rev. Lett. 

46, 614 (1981). 

31 K. Machida, J. Phys. Soc. Jpn. 50, 2195 (1981). 

32 A. M. Gabovich and A. S. Shpigel, J. Phys. F: Met. Phys. 
14, 3031 (1984). 

33 M. Gulacsi and Zs. Gulacsi, Phys. Rev. B 33, 6147 (1986). 
3* M. Kato and K. Machida, Phys. Rev. B 37, 1510 (1988). 
3^ M. L. Kulic, E. Goreatchkovski, A. I. Lichtenstein, and M. 

Mehring, Physica C 252, 27 (1995). 
3® W. A. Atkinson, Phys. Rev. B 75, 024510 (2007). 
3^" T. Das, R. S. Markiewicz, and A. Bansil, Phys. Rev. B 77, 



134516 (2008). 
3^ A. B. Vorontsov, M. G. Vavilov, and A. V. Chubukov, 

Phys. Rev. B 79, 060508(R) (2009). 
3^ V. Stanev, J. Kang, and Z. Tesanovic, Phys. Rev. B 78, 

184509 (2008). 
^° D. Parker, M. G. Vavilov, A. V. Chubukov, and I. I. Mazin, 

Phys. Rev. B 80, 100508(R) (2009). 
"1 M. G. Vavilov, A. V. Chubukov, and A. B. Vorontsov, 

ar Xiv:09 12.3556 (2009). 
"2 A. B.Vo rontsov, M. G. Vavilov, and A. V. Chubukov, 

arXiv:l003.2389 (2010). 

"3 P. Ghaemi and A. Vishwanath. larXiv:1002.4638l 

"" L I. Mazin, D. J. Singh, M. D. Johannes and M. H. Du, 

Phys. Rev. Lett. 101, 057003 (2008). 
"5 K. Kuroki, S. Onari, R. Arita, H. Usui, Y. Tanaka, H. Kon- 

tani, and H. Aoki , Phys. Rev. Lett. 101, 087004 (2008). 
*® A.V. Chubukov, D. V. Efremov, and L Eremin, Phys. Rev. 

B 78, 134512 (2008). 
"■^ H. Ikeda, J. Phys. Soc. Jpn. 77, 123707 (2008). 
"^^ F. Wang, H. Zhai, Y. Ran, A. Vishwanath, and D.-H. Lee, 

Phys. Rev. Lett. 102, 047005 (2009). 
*® R. Sknepnek, G. Samolyuk, Y. B. Lee, and J. Schmalian, 

Phys. Rev. B 79, 054511 (2009). 
5° S. eraser, T. A. Maier, P. J. Hirschfeld, and D. J. 

Scalapino, New J. Phys. 11, 025016 (2009). 
^1 J. Zhang, R. Sknepnek, R. M. Fernandes, and J. 

Schmahan, Phys. Rev. B 79, 220502(R) (2009). 
^2 L L Mazin and J. Schmalian, Physica C 469, 614 (2009). 
^3 C. Liu,T. Kondo, R. M. Fernandes, A. D. Palczewski, E. D. 

Mun, N. Ni, A. N. Thaler, A. Bostwick, E. Rotenberg, J. 

Schmalian, S. L. Bud'ko, P. C. Canfield, and A. Kaminski, 

arXiv:0910.1799 (2009). 
^^ D. J. Van Harlingen, Rev. Mod. Phys. 67, 515 (1995). 
^^ C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 

(2000). 
^^ C. T. Chen, C. C. Tsuei, M. B. Ketchen, Z. A. Ren, Z. X. 

and Zhao, Nature Phys. 6, 260 (2010). 
^'^ Y. Sekiba, T. Sato, K. Nakayama, K. Terashima, P. 

Richard, J. H. Bowen, H. Ding, Y.-M. Xu, L. J. Li, G. 

H. Cao, Z.-A. Xu, and T. Takahashi, New J. Phys. 11, 

025020 (2009). 
^* V. Brouet, M. Marsi, B. Mansart, A. Nicolaou, A. Taleb- 

Ibrahimi, P. Le Fevre, F. Bertran, F. Rullier-Albenque, A. 

Forget, and D. Colson, Phys. Rev. B 80, 165115 (2009). 
'^^ C. Liu, A. D. Palczewski, T. Kondo, A. N. Thaler, N. Ni, S. 

L. Bud'ko, P. C. Canfield, and A. Kaminski, unpublished. 
®° Notice that this quantity is a factor of 2 greater than the 

g we defined in our previous work, Ref. 20. 
^^ J. M. Kosterlitz, D. R. Nelson, and M. E. Fisher, Phys. 

Rev. B 13, 412 (1976). 
^2 A. Aharony, J. Stat. Phys. 110, 659 (2003). 
"^3 A. Kreyssig, M. G. Kim, S. Nandi, D. K. Pratt, W. Tian, J. 

L. Zarestky, N. Ni, A. Thaler, S. L. Bud'ko, P. C. Canfield, 

R. J. McQueeney, and A. L Goldman, Phys. Rev. B 81, 

134512 (2010). 
^^ S.-C. Zhang, Science 275, 1089 (1997). 
^^ E. Demler, W. Hanke, and S.-C. Zhang, Rev. Mod. Phys. 

76, 909 (2004). 
^^ D. Podolsky, H.-Y. Kee, and Y. B. Kim, Europhys. Lett. 

88, 17004 (2009). 
^■^ L. V. Keldysh and Y. V. Kopaev, Sov. Phys. Solid State 

6, 2219 (1965). 
^^ J. des Cloizeaux, J. Phys. Chem. Solids 26, 259 (1965). 
^3 P. C. Canfield and G. W. Crabtree, Phys. Today 56, 34 



24 



(2003). 
™ P. Hohenberg, Phys. Rev. 158, 383 (1967); N. D. Mermin 

and H. Wagner, Phys. Rev. Lett. 17, 1133 (1966). 
■^^ D. Parker and I. I. Mazin, Phys. Rev. Lett. 102, 227007 

(2009). 
^2 J. Wu and P. Phillips, Phys. Rev. B 79, 092502 (2009). 
'^^ J. Linder, L Sperstad, and A. Sudb0, Phys. Rev. B 80, 

020503(R) (2009). 
■^^ A. V. Chubukov, Physica C 469, 640 (2009). 
'^^ A. V. Chubukov, M. G. Vavilov, and A. B. Vorontsov, 

Phys. Rev. B 80, 140515(R) (2009). 
''^ R. Thomale, C. Piatt, J. Hu, C. Honerkamp, and B. A. 

Bernevig, Phys. Rev. B 80, 180505(R) (2009). 
"^"^ P. Bonville, F. Rullier-Albenque, D. Colson, and A. Forget, 

larXiv:1 002.0931 (2010). 
■^® T. M. Rice, Phys. Rev. B 2, 3619 (1970). 
"^^ Y. Zhang, E. Dernier, and S. Sachdev, Phys. Rev. B 66, 

094501 (2002). 
®° C. Fang, H. Yao, W.-F. Tsai, J. Hu, and S. A. Kivelson, 

Phys. Rev. B 77, 224509 (2008). 
*i C. Xu, M. Miiller, and S. Sachdev, Phys. Rev. B 78, 

020501 (R) (2008). 
*^ R. M. Fernandes, L. H. VanBebber, S. Bhattacharya, P. 

Chandra, V. Keppens, D. Mandrus, M. A. McGuire, B. 

C. Sales, A. S. Sefat, and J. Schmahan, ,arXiv:0911. 30841 

(2009). 
*^ For a recent review, see G. Knebel, D. Aoki, and J. Flou- 

quet, arXiv:0911.5223 (2009). 
*" E. D. Isaacs, P. Zschack, C. L. Broholm, C. Burns, G. 

Aepph, A. P. Ramirez, T. T. M. Palstra, R. W. Erwin, N. 



Stiicheli, and E. Bucher, Phys. Rev. Lett. 75, 1178 (1995). 
**^ S. Sachdev, arXiv: 1002.3823 (2010). 
*** S. Sachdev, Quantum phase transitions, Cambridge Univ. 

Press (2001). 
*^ P. C. Canfield, P. L. Gammel, and D. J. Bishop, Phys. 

Today 51, 40 (1998). 

For a review, see Superconductivity in Ternary Compounds, 

Vol. II, edited by M. B. Maple and 0. Fisher (Springer 

Verlag, Berlin, 1982). 
^^ P. Coleman, in Handbook of Magnetism and Advanced 

Magnetic Materials, edited by H. Kronmuller and S. 

Parkin, Vol. 1, 95 (Wiley, New York, 2007). 
^° M. Rotter, M. Tegel, L Schellenberg, F. M. Schappacher, 

R. Pottgen, J. Deisenhofer, A. Giinther, F. Schrettle, A. 

Loidl, and D. Johrendt, New J. Phys. 11, 025014 (2009). 
^^ Y. Kamihara, T. Nomura, M. Hirano, J. E. Kim, K. Kato, 

M. Takata, Y. Kobayashi, S. Kitao, S. Higashitaniguchi, Y. 

Yoda, M. Seto, and H. Hosono, New J. Phys. 12, 033005 

(2010). 
^2 R. T. Gordon, N. Ni, C. Martin, M. A. Tanatar, M. D. Van- 

nette, H. Kim, G. D. Samolyuk, J. Schmalian, S. Nandi, 

A. Kreyssig, A. L Goldman, J. Q. Yan, S. L. Bud'ko, P. C. 

Canfield, and R. Prozorov, Phys. Rev. Lett. 102, 127004 

(2009). 
^^ S. Nandi, M. G. Kim, A. Kreyssig, R. M. Fernandes, D. 

K. Pratt, A. Thaler, N. Ni, S. L. Bud'ko, P. C. Canfield, 

J. Schmalian, R. J. McQueeney, and A. L Goldman, Phys. 

Rev. Lett. 104, 057006 (2010). 



